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GENERAL  RELATIVITY 
AND  ITS  EFFECTS  ON  PLANETARY  ORBITS 
AND  INTERPLANETARY  OBSERVATIONS 


I.  INTRODUCTION 

New  radar  installations,  such  as  the  7840-Mcps  Haystack  facility  developed  by  the  M.I.T. 
Lincoln  Laboratory,  will  shortly  provide  much  more  data  on  the  time  delays  and  Doppler  shifts 
of  interplanetary  radar  signals.  These  data,  when  combined  with  other  optical  and  radar  obser¬ 
vations,  will  allow  some  of  the  predictions  of  Einstein’s  General  Relativity  Theory  (GRT)  to  be 
tested.  Conventional  astronomy  assumes  that  the  equations  governing  the  motion  of  celestial 
bodies  are  those  of  Newtonian  Gravitational  Theory  (NGT),  and  that  those  governing  the  propa¬ 
gation  of  electromagnetic  signals  are  given  by  the  Special  Relativity  Theory  (SRT).  The  corre¬ 
sponding  equations  derived  from  GRT  are  conventional  in  the  lowest  approximation,  but  small 
correction  terms  appear  in  higher  approximations;  the  tests  will  be  performed  to  determine 
the  presence  or  absence  of  these  terms. 

In  this  report,  we  first  discuss  the  assumptions  that  underlie  GRT  and  introduce  the  mathe¬ 
matical  concepts  and  techniques  that  are  to  be  used.  We  then  obtain  the  approximate  equations 
of  motion  to  the  accuracy  needed  for  the  solar  system.  Since  there  is  a  good  deal  of  contro¬ 
versy  among  workers  in  the  field  with  regard  to  how  the  approximate  equations  are  to  be  obtained 
from  the  exact  ones,  the  views  and  assumptions  of  the  various  schools  are  discussed  at  some 
length.  All  the  approaches  yield  the  same  equations  of  motion  for  the  solar  system  except  for 
higher-order  terms  which  are  so  small  as  to  be  unmeasurable.  The  propagation  of  electro¬ 
magnetic  waves  according  to  GRT  is  then  considered,  and  expressions  for  the  time-delay  and 
Doppler-shift  corrections  obtained. 
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II.  THE  FIELD  EQUATIONS  AND  EQUATIONS  OF  MOTION 


Let  us  begin  by  writing  down  the  field  equations  of  CRT  and  discussing  the  expressions  in 
them: 1 


r..  -  j  g-u  R  =  t  l, 

ik  2  foik  4  ik 

c 


(2-1) 


In  explanation  of  the  notation,  it  must  be  understood  that  we  are  dealing  with  functions  in  a 
four-dimensional  space  which  has  the  property  of  being  locally  Minkowskian,  i.e.,  we  can,  by 
using  suitably  chosen  coordinates,  make  an  infinitesimally  small  neighborhood  around  each  point 
look  like  the  usual  space-time  of  special  relativity.  If  we  have  an  expression  written  down  in  one 
coordinate  system,  then  the  indices  tell  us  how  to  rewrite  it  in  another  coordinate  system,  the 
general  rule  being^ 
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where  xl  and  x1  are  the  coordinates  of  the  same  physical  point  in  the  two  systems.  In  Eq.  (2-2) 
and  in  the  sequel  we  use  the  Einstein  summation  convention:  A  repeated  latin  index  is  to  be 
summed  from  one  to  four.  Objects  transforming  according  to  Eq.  (2-2)  are  called  tensors  and 
because  of  the  linearity  of  the  transformation  law,  a  tensor  equation  has  the  same  form  in  all 
coordinate  systems.  Thus  Eq.  (2-1)  is  shorthand  for  a  set  of  16  equations  and  can  be  written 
down  in  any  coordinate  system. 

SRT  rests  essentially  on  three  postulates: 


(a)  There  is  a  class  of  privileged  reference  frames,  called  inertial  frames, 
in  which  the  laws  of  physics  are  the  same. 

(b)  The  velocity  of  a  light  signal  appears  to  be  the  same  when  measured 
in  each  inertial  frame. 

(c)  Space-time  is  homogeneous. 

Translating  this  into  mathematical  language,  we  find  that  the  Cartesian  coordinates  assigned  to 
an  event  in  space-time  in  different  inertial  frames  are  connected  by  Lorentz  transformations, 
and,  in  order  to  have  the  equations  retain  their  form,  that  basic  equations  should  be  tensor 
equations  (in  SR  we  do  not  deal  with  general  nonsingular  transformations  of  the  form  of 
Eq.  (2-3)  but  only  with  Lorentz  transformations).  A  Lorentz  transformation  is  a  linear  trans¬ 
formation  with  the  property  that 
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if  A  and  B  are  two  events.  By  looking  at  two  infinitesimally  close  events  we  can  rewrite 
Eq. (2-4)  as 
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Therefore,  Lorentz  transformations  are  linear  transformations  which  leave  the  quadratic  form 
in  Eq.  (2-5)  invariant. 

We  can  also  formulate  Lorentz  transformations  in  a  general  curvilinear  coordinate  system, 
provided  we  know  the  relationship  between  the  general  coordinates  and  a  set  of  inertial  Cartesian 
coordinates,  by  writing 
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(X  are  inertial  Cartesian  coordinates) 


in  which  case 
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and  all  tensor  equations  can  likewise  be  modified  trivially.  Finally,  we  note  that,  relative  to  an 

inertial  system,  a  free  particle  moves  in  a  straight  line  in  space-time.  Since  a  straight  line 

can,  in  ordinary  space,  be  characterized  by  a  variational  principle  (the  arc  length  between  two 

points  is  shortest  along  the  straight  line  connecting  them),  we  might  expect  a  similar  result  to 
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hold  in  Minkowski  space.  It  turns  out  that  if  we  demand 


rB 

6  \  ds  =  0  (2-11) 
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where  A  and  B  are  two  points  on  the  particle's  trajectory,  and 

ds2=-ij..dXldXj  (2-12) 

and  use  the  standard  techniques  of  the  calculus  of  variations,  then  we  obtain  as  Euler- Lagrange 
equations 

d2Xl  _  n  (2-13) 
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ds 


If  we  transform  to  curvilinear  coordinates  and  use  Eq.(2-9)  and  the  variational  principle,  then 
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in  place  of  Eq.  (2-13)  we  obtain 
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glmgmk  =  6k  ’  6k  =  diagonal  (1,  1,  1,  1)  .  (2-16) 

Note  that  if  g^  =  77^  then  Eq.  (2-14)  reduces  to  Eq.  (2-13).  A  trajectory  which  satisfies  Eq.  (2-14) 

is  called  a  geodesic.  The  are  called  Christoffel  symbols  of  the  second  kind;  they  are  not 

quite  tensors  (if  they  were,  they  would  vanish  in  all  coordinate  systems  if  they  vanished  in  one) 

but  their  transformation  properties  are  of  no  interest  to  us.  If  we  want  the  path  of  a  light  signal, 

we  cannot  use  ds  as  a  parameter  since  ds  =  0  along  the  trajectory;  it  can  be  shown  that  Eq.  (2-14) 
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can  then  be  replaced  by 


d2x1  r  i  dxk  dx*  _  dxj_  dxk 

dx2  kf  d\  d\  ;  gik  dx  dx 


(2-17) 


Equations  (2-17)  define  a  null  geodesic;  they  determine  the  path  and  the  parameter  (called  a 
"distinguished  parameter  along  the  geodesic"). 

GRT  holds  that  if  matter  is  present,  particles  still  move  along  geodesics,  but  the  g.  no 
longer  have  the  form  of  Eq.  (2-10)  so  the  geodesics  are  not  straight  lines.  As  stated  above,  we 
assume  space-time  is  locally  Minkowskian  so  that  it  is  possible  to  introduce  coordinates  in 
terms  of  which 
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where  P  is  any  point,  but  in  these  coordinates  g.^(Q)  =£  77  in  general  if  Q  is  not  in  an  infinites¬ 
imal  neighborhood  of  P.  If  it  is^  possible  to  introduce  a  coordinate  system  in  which  g. 
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globally  the  space  is  called  flat.  As  we  will  see,  the  left-hand  side  of  the  field  equations  (2-1) 
vanishes  for  a  flat  space. 

Given  a  set  of  sixteen  functions  g^(x)  with  g-  =  g^  how  can  we  tell  if  the  space  is  flat?  This 
is  a  problem  in  partial  differential  equations,  since  it  is  obviously  equivalent  to  asking  under 
what  circumstances  the  equations 
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[Eq.  (2-10)] 


have  solutions  for  X1.  The  problem  is  quite  analogous  to  the  question  as  to  whether  a  given 
three-dimensional  vector  field  can  be  written  as  the  gradient  of  a  scalar: 
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The  necessary  and  sufficient  condition  for  Eq.(2-19)  to  be  true  is  known  from  elementary 
vector  analysis  to  be 
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The  corresponding  necessary  and  sufficient  condition  for  Eq.(2-10)  to  be  soluble  is^’^ 
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Rikim  is  called  the  Ri^mann  or  curvature  tensor.  The  term  "flat  space"  and  "curvature"  come 
from  the  study  of  two-dimensional  surfaces  in  ordinary  three-dimensional  space.  For  coordi¬ 
nates  on  a  curved  surface,  the  two-dimensional  Riemann  tensor  is  proportional  to  the  ordinary 
curvature  (Gaussian  curvature)  of  the  surface  7  The  tensors  appearing  on  the  left-hand  side  of  the 

g 

field  equations  (2-1)  are  defined  as 
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Because  of  the  way  it  is  defined,  the  Riemann  tensor  has  a  great  deal  of  symmetry  (it  has  only 
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20  independent  components  instead  of  4  =  256)  and  its  derivatives  also  satisfy  some  identities 

(Bianchi  identities).  The  situation  is  again  analogous  to  the  case  of  the  curl  of  a  vector  in  three 
dimensions  where,  from  Eq.  (2-20)  it  can  be  seen  that 
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and 


V  .  (V  X  A)  =  0  (2-26) 

with  Eq.  (2-26)  being  identically  true,  whether  or  not  V  x  A  vanishes.  If  we  define 
def  i 

G.^  =  R^  —  (Einstein  tensor)  (2-27) 

and 
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then  the  Bianchi  identities  take  the  form 
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It  is  very  important  to  notice  that  Eq.  (2-29)  is  identically  true  because  of  the  way  the  Einstein 
tensor  is  related  to  the  Riemann  tensor.  We  accordingly  conclude  that  a  necessary  condition 
for  the  field  equations  of  GRT  to  be  consistent  is 
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Note:  It  is  conventional  to  consider  two  tensors  which  are  related 
as  in  Eq.  (2-31)  to  be  the  "same"  tenspr  and  to  speak  of  lowering  and 
raising  indices  by  means  of  g. .  and  g1!,  thus,  e.g., 
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ij  &ik  j  6ik6ji 

The  same  basic  symbol  is  used  for  tensors  related  in  this  manner. 

T.^  is  a  generalization  of  the  energy-momentum  tensor  of  special  relativity.  In  flat  space, 
if  we  have  a  tensor  which  satisfies  the  condition 
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where  x  are  Cartesian  coordinates,  then 
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and  integrating  over  all  space  between  two  times,  we  find 


(2-33) 
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ik 

so  if  T  vanishes  sufficiently  rapidly  at  spatial  infinity,  then  the  four  quantities  P.  are  con- 

1  ik 

served.  If  a  theory  is  based  on  a  variational  principle  and  is  Lorentz-invariant,  such  a  T  can 
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always  be  found  if  no  "external  fields"  are  present.  This  situation  is  reminiscent  of  the  re¬ 
sult  of  nonrelativistic  particle  dynamics,  where  the  energy  and  momentum  of  a  system  of  parti¬ 
cles  are  conserved  if  no  external  forces  act  on  it;  it  can  be  shown  that  it  is  consistent  to  inter- 

i  ik 

pret  P  as  the  energy-momentum  vector  of  the  relativistic  field  (the  tensorial  nature  of  T  will 

guarantee  that  the  P.  transform  like  a  four-vector  under  Lorentz-transformation),  i.e.,  P  is 

the  momentum  and  E  the  energy. 

Suppose  we  wish  to  discuss  "particles."  A  "particle,"  by  definition,  has  a  location,  no 
extension  in  space  and  is  completely  characterized  by  a  parameter  " m ^ "  called  its  rest  mass, 
and  by  its  "velocity."  It  has  an  energy  and  momentum  related  by 

t-,  2  2  ^  224 

E  -  c  P  =  m  c  .  (2-37) 

o 


6 


1 1 

We  are  thus  more  or  less  forced  to  take  for  T..  the  form 

lk 

Tlk  =  m  c  —■  6(r  -  R)  ,  ds2  =  ( 1  -  v2/c2)  dt2  (SRT)  (2-38) 

o  as  dt 

This  leads  to 

=  (m  yv.m  cy)  ( y  = - - -  |  (2-39) 

°  °  V 

which  is  the  familiar  result  of  relativistic  mechanics. 

How  seriously  should  one  take  the  idea  of  "point"  particles?  In  electrodynamics,  the  use 
of  6-functions  leads  to  infinite  self-energies  and  other  calamities  which  have  to  be  argued  away 
by  "renormalization"1^  and  which  make  the  theory  formally  inconsistent  (this  formal  inconsist¬ 
ency  may  just  be  due  to  not  using  a  sophisticated  enough  mathematical  formalism  for  handling 

6-like  functions,13  but  no  one  really  knows).  On  the  other  hand,  for  extended  matter  distribu- 

ik 

tions,  no  simple  form  for  T  is  possible  since  extended  matter  must  involve  internal  degrees 

of  freedom  ("rigid"  bodies  can  only  exist  if  there  are  forces  which  act  instantaneously  at  a  dis- 

1 4 

tance  and  this  violates  the  postulates  of  SRT).  People  who  are  unhappy  about  using  6-functions 
but  do  not  want  to  include  internal  forces  explicitly,  usually  just  let 


m  6(7’  —  R)  p(7)  (2-40) 

o 

-►  15 

where  p(r)  is  "concentrated"  around  R,  and  take  the  6-function  limit  wherever  it  makes  sense. 
This  is  called  the  energy-momentum  tensor  for  "incoherent  matter."  In  order  to  include  the 
interactions  between  the  different  parts  of  an  extended  body,  elasticity  theory  and  thermo¬ 
dynamics  must  be  used  as  a  guide  to  determine  the  phenomenological  form.1  17  It  is  not  clear 
how  a  more  fundamental  level  than  the  macroscopic  one  can  be  reached  with  this  approach,  but 
on  the  other  hand  it  may  very  well  be  that  all  classical  theories  are  inherently  macroscopic  since 
quantum  theory  becomes  important  for  elementary  particles  and  for  very  small  distances.  For 
astronomical  problems  the  various  bodies  are  so  far  away  that  internal  forces  can  be  effectively 
ignored  and  the  forces  between  celestial  objects  can  be  calculated  as  if  there  were  point  planets 
(provided  "renormalization"  takes  place,  since  for  self-interactions  the  particle  approximation 
is  invalid);  in  the  language  of  electrodynamics,  a  multipole  expansion  is  used,  keeping  only  the 

monopole  moment  (the  problem  of  incorporating  higher  moments,  e.g.,  spinning  particles,  has 

1 8 

not  yet  been  solved  on  a  nonphenomenological  level). 

19 

We  now  come  to  one  of  the  most  remarkable  results  in  general  relativity.  In  1927  Einstein 
and  Grommer  found  that  the  equations  of  motion  of  an  isolated  test  particle  could  be  derived 
from  the  field  equations.  To  obtain  this  we  first  rewrite  Eq.(2-40)  to  make  it  a  tensor  under  a 
general  transformation;  the  change  is  trivial:^ 
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We  then  integrate  Eq.  (2-30)  around  a  tube  in  space-time  which  encloses  the  world-line  of  the 
test  particle  and  no  other  matter,  although  other  matter  may  be  present  outside  of  the  tube. 

The  properties  of  the  6-function  under  integration  then  lead  to  the  result  that  the  particle  obeys 
the  geodetic  equations  of  motion: 


,2  i  .  .  i  .  k 

d  x  r  1  dx  dx 

.2  ki  ds  ds 

ds 


(2-43) 


There  has  been  one  bit  of  ’’swindling"  though;  T  as  defined  in  Eq.  (2-41)  depends  explicitly 
on  g.j,  and  if  we  have  a  particle  we  would  expect  it  to  produce  a  metric  tensor  which  is  singular 
at  the  particle  (analogous  to  the  Coulomb  potential  at  a  point  charge).  The  product  of  a  singular 
function  and  a  6-function  must  be  defined  more  carefully  and  in  fact  may  be  impossible  to  define. 
It  is  for  this  reason  that  we  call  the  particle  a  "test”  particle  and  assume  that  it  does  not  con¬ 
tribute  to  the  field  so  that  the  coefficient  of  6(r—  R)  in  Eq.(2-41)  is  well  defined.  "Test"  par¬ 
ticles  make  intuitive  sense  but  may  not  make  mathematical  sense  in  a  nonlinear  theory  like 
GRT  where  the  superposition  principle  does  not  hold.  Nevertheless  one  is  usually  forced  to 

think  in  terms  of  test  particles  in  order  to  relate  the  theory  to  pre-relativity  physics. 

2 1 

How  can  we  obtain  Newtonian  Gravitational  Theory  (NGT)  from  GRT?  We  would  expect 
NGT  to  apply  when  we  can  ignore  the  finite  velocity  of  light  both  for  test  particles  and  for  the 
sources  of  the  field  as  well  as  the  nonlinearity  of  the  field  ("weak  field")  since  NGT  is  a  linear 
theory.  We  therefore  assume  that  it  is  possible  to  find  a  coordinate  system  in  which  the  metric 
can  be  expanded  into  a  series 


gu" 


ij 


+  \h. . 
ilJ 


+  \2h. . 


(2-44) 


where  X  is  a  formal  parameter,  and  we  keep  terms  to  order  X.  We  further  assume  that  in 
this  coordinate  system  the  time  dependence  of  g^.  is  entirely  due  to  the  motion  of  the  sources. 
This  corresponds,  intuitively,  to  assuming  there  is  no  "gravitational  radiation"  incident  on  the 
system  which  is  therefore  "isolated."  (The  term  "gravitational  radiation"  is  in  quotes  because 
there  is  a  great  deal  of  controversy  over  what  is  meant  by  radiation  in  the  theory  and  what 
effect,  if  any,  it  would  have  on  the  motion  of  a  system. )^  Since  the  velocities  of  the  sources 
are  assumed  small 
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where  v  is  some  characteristic  velocity  associated  with  the  matter.  For  the  same  reason  the 
equations  of  motion  of  a  test  charge  can  be  simplified  by  using 


4 

dx 

~  «  1  = 
c 

dx 

cdt 

cdt 
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The  equations  of  motion  become,  approximately 


(2-47) 
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(2-48) 
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which  is  to  be  compared  to  Newton's  Law 


M-=-V(p  .  (2-49) 

dt 


In  order  to  find  the  relation  between  hAA  and  the  sources  of  the  field,  we  use  the  T.,  due  to 

1 44  lk 

incoherent  matter  and  the  small  velocity  approximation, 

T44  «p  »  (other  components)  (2-50) 


where  p  is  the  density  of  matter,  and  obtain24 


VZ^44--HG»  ■ 

1  c 


We  also  identify 


g44 


a  -1  - 


2  (p 

“Z 
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V  (p  =  +4t rGp 


(2-51) 


(2-52) 

(2-53) 


which  is  just  the  equation  for  the  Newtonian  potential  if  G  is  the  usual  gravitational  constant. 
The  form  of  the  results  is  most  gratifying  but  it  is  very  important  to  examine  closely  the  extent 
to  which  we  have  derived  Newtonian  Mechanics.  We  therefore  note  the  following  points: 


(a)  We  have  not  solved  the  field  equations  but,  rather,  the  formally 

linearized  equations  obtained  by  writing  g^  =  77^  +  \h.  ^  and  dropping 

terms  quadratic  or  higher  in  X.  It  is  not  known  to  what  extent  there 
is  a  correspondence  between  solutions  of  the  linearized  equations 
and  of  the  full  equations.  Solutions  of  the  linearized  equations  can 
obviously  be  superposed  to  give  a  new  solution,  which  is  not  generally 
true  of  the  nonlinear  full  equations;  on  the  other  hand,  the  full  equa¬ 
tions  are  covariant  under  arbitrary  coordinate  transformations  and 
the  linearized  ones  are  not.^5 


(b)  NGT  applies  to  bodies  moving  at  arbitrary  velocities  in  arbitrarily 
strong  gravitational  fields,  and  the  basic  equations  (2-53)  and  (2-49) 
reflect  this.  There  is  nothing  in  these  equations  themselves  that 

forces  us  to  use  solutions  with  v/c  small  or  for  which  |2<p/c2|«  1. 
Thus  there  are  solutions  of  the  approximate  equations  which  violate 
the  conditions  required  to  derive  the  approximate  equations  from  the 
full  theory. 

(c)  For  the  N-body  system  the  Newtonian  potential,  ignoring  structure,  is 


N 

<p(r)  =  -G  Yj 


This  serves  to  define  the  "mass"  of  the  ith  body  in  NGT  which  must, 
of  course,  be  determined  experimentally  by  fitting  parameters  to  a 
solution,  as  must  the  gravitational  constant.  In  the  first  approxima¬ 
tion  to  GRT  we  end  up  with  the  same  potential  by  identifying  with 
the  integral  over  the  body  of  the  p  in  Eq.  (2-53).  In  higher  approxi¬ 
mations,  however,  there  will  be,  among  other  changes,  corrections 
to  cp  which  will  make  it  nonlinear  in  the  individual  masses.  Since 


(2-54) 
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the  "Newtonian"  mass  is  determined  as  that  value  of  a  parameter 
which  gives  a  best  fit  to  a  solution  with  (p  in  the  form  of  Eq.  (2-54), 
we  may  expect  the  masses  determined  from  GRT  to  differ  from 
the  "Newtonian"  values. 


Before  discussing  higher  approximations  to  the  field  equations,  it  is  convenient  to  have  on 

hand  a  particular  exact  solution  of  the  field  equations,  the  Schwarzschild  solution.  This  solu- 

..  .  26 
tion  is 


ds  = 


2  ,+2 
c  dt  — 


?  ?  ?  2  2  2 
(dr  +  r  d0  +  r  sin  0d<p  ) 


(2-55) 


where 


2Gm 


2 

c 


(2-56) 


What  physical  situation  does  this  correspond  to?  If  we  substitute  the  g^  of  Eq.  (2-55)  into 
the  field  equations  (2-1)  we  find  vanishes  everywhere  except  at  the  point  r  =  0  when  the  g.^ 
and  their  derivatives  are  singular.  Thus  the  solution  represents  matter  at  the  origin.  In  fact, 
this  solution  could  be  obtained  if  the  most  general  modification  of  the  Minkowski  metric  in 
spherical  coordinates  were  demanded,  i.e., 

222  22222  2 

ds  =  c  dt  —(dr  +  r  d0  +  r  sin  0d <p  )  (2-57) 

which  would  not  destroy  its  spherical  symmetry  in  space  and  its  symmetry  with  relation  to  past 
27 

and  future  time.  The  general  solution,  up  to  a  coordinate  transformation,  is  that  of  Eq.  (2-55). 

The  quantity  "r  "  appears  as  a  constant  of  integration  and  its  identification  with  the  "mass"  of 

°  2 
the  matter  at  the  origin  is  made  by  examining  the  expression  for  ds  when  r  »  rQ, 

ds2  «  (l  -fr)  c2dt2-(l  +  (dr2  +  r2de2  +  r2  sin20d«>2)  (2-58) 

and  appealing  to  the  correspondence  with  the  Newtonian  potential  as  given  by  Eq.  (2-52).  It  is 
therefore  said  that  the  Schwarzschild  solution  represents  the  field  of  a  structureless  point  mass 
fixed  at  the  origin. 

One  must  not  be  misled  by  the  fact  that  we  use  the  same  symbols  for  the  coordinates, 
viz.  (r,  0,  cp ,  t),  in  the  Minkowski  metric,  in  the  Schwarzschild  metric,  and  in  the  Newtonian 
equations  of  motion.  The  Schwarzschild  metric  is  that  of  a  non-Euclidean  space  and  there  is 
no  quantity  which  has  all  the  properties  of  the  Euclidean  radius  vector  (i.e.,  being  equal  both 
to  the  distance  from  the  center  and  to  the  length  of  the  circumference  of  a  circle  divided  by 
2 7r ) .  By  changing  the  coordinates  in  which  we  express  the  Schwarzschild  solution,  we  can  pro¬ 
duce  an  infinite  number  of  different  "radius  vectors,"  each  of  which  can  be  identified  with  the 
Newtonian  one  far  away  from  the  source;  all  we  must  guarantee  is  that  the  asymptotic  g. .  be  of 
the  form 

g .. - 1-—  ,  g  0  ~  >5  „  +  0  (— )  (2-59) 

fc*44  r  »  a  p  r 

and  that  the  relationship  between  the  old  and  the  new  variables  be  in  some  sense  "smooth."  The 
whole  question  of  which  coordinates  to  use  to  express  the  Schwarzschild  solution  is  a  thorny  one. 
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The  coordinates  used  in  Eq.  (2-55)  are  the  so-called  "isotropic"  ones^  and  it  will  be  noted  that 
g^  =  0  for  r  =  r  /4,  the  so-called  " Schwarzschild  radius."  There  is  a  strong  temptation  to  ig¬ 
nore  this  since  for  all  macroscopic  bodies  rQ  lies  well  in  the  interior  where  it  would  not  be 
expected  that  the  Schwarzschild  solution  would  be  a  good  approximation  to  the  real  solution;  it 
is,  however,  at  least  esthetically  disturbing.  It  has  been  found  that  the  peculiar  behavior  at 
r  /4  is  due  to  the  choice  of  coordinates  and  can  be  eliminated  by  using  different  ones,  but  then 

other  strange  things  happen  to  light  signals  which  penetrate  the  barrier.  How  much  of  this  is 

28 

mathematics  and  how  much  is  physics  is  not  known  as  yet. 

The  form  of  the  equations  of  the  geodesics  will  be  different  for  each  choice  of  "r,"  but,  pro¬ 
vided  that  Eq.  (2-59)  holds  and  that  no  spurious  zeros  or  infinities  are  introduced  by  the  change 

2 

of  variables,  the  force  will  look  like  a  "Newtonian"  l/r  one  plus  additional  terms  which  will  be 

29 

small  if  "r"  is  identified  with  the  "Newtonian"  radial  variable.  As  Eddington  points  out,  GRT 
makes  perfectly  definite  predictions;  the  ambiguity  lies  in  Newtonian  theory  which  assumes  that 
the  force  is  strictly  l/r  in  nature.  In  GRT,  "r"  is  just  a  marker.  If  we  just  use  the  asymptotic 
form  (taking  the  effects  of  planet-planet  interactions  into  account  in  a  manner  that  amounts  to 
using  Newton’s  Law  according  to  Eq.  (2-52),  then  we  have  exactly  the  usual  Newtonian  mechanics. 
The  extra  terms  for  the  sun's  field  produce  the  famous  advance  of  the  perihelia  of  the  planets. 

The  precession  of  an  almost  closed  orbit  clearly  has  an  observable  meaning  which  is  relatively 
independent  of  the  choice  of  "r, "  so  different  "r"s  give  the  same  approximate  formulas.  This 
mixture  of  GRT  and  NGT  is,  as  is  well  known,  in  very  good  agreement  with  observation. 

This  ambiguity  in  "r,"  about  which  much  has  been  written  (much  of  it  wrong!)  is  in  fact  to 

25 

be  expected  and  should  not  bother  us  unduly.  As  Bergmann  points  out,  the  fact  that  the  equa-- 
tions  of  GRT  look  the  same  in  any  coordinate  system  is  not  just  a  minor  generalization  on  equa¬ 
tions  looking  the  same  in  any  Galilean  or  Minkowskian  inertial  frame.  Mathematically  the  dif¬ 
ference  is  between  invariance  under  an  infinite  dimensional  Lie  group  and  invariance  under  a 
finite  dimensional  Lie  group.  The  analogue  of  general  invariance  is  gauge  invariance  in  elec¬ 
tromagnetic  theory,  the  coordinates  in  GRT  being  the  analogues  of  the  potentials  in  electro¬ 
dynamics.  Given  the  electrodynamic  potentials,  we  certainly  can  calculate  the  fields,  but  since 
all  we  can  observe  in  nature  are  gauge- invariant  quantities,  we  never  can  determine  the  poten¬ 
tials  experimentally. 

Much  work  has  gone  into  the  problem  of  determining  a  "complete  set  of  observables"  for 

25 

GRT,  but  little  has  been  accomplished.  No  one  has  observed  "r"  for  a  celestial  body;  w’hat 
has  been  observed  is  light  or  radar  reflected  from  a  planet  or,  say,  emitted  by  a  comet.  GRT 
tells  us  how  this  electromagnetic  radiation  propagates  and  how  it  is  affected  by  its  passage 
through  a  gravitational  field.  NGT  makes  no  predictions  along  these  lines.  For  different  choices 
of  "r,"  the  behavior  of  light  rays  will  be  different  but  quantities  like  time  delays  have  an  unambig¬ 
uous  meaning  since  they  measure  the  time  it  takes  for  a  signal  to  travel  from  one  point  on  the 
earth  to  the  planet  and  back  again.  In  different  coordinate  systems  we  would  assign  different 
values  of  the  "markers"  to  the  positions  of  the  earth  and  of  the  planet  relative  to  the  sun,  but 
the  relationship  between  the  assigned  "markers"  and  the  time  delay  is  unambiguous  and  can  be 
checked  experimentally. 

We  shall  discuss  optics  in  the  Schwarzschild  field  further  in  a  later  section  of  this  report. 

Let  us  now  turn  to  the  problem  of  the  equations  of  motion  for  an  N-body  system.  If  the  solar 
system  is  considered  isolated,  then  the  usual  Newtonian  equations  determine  the  motion  of  the 
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planets  relative  to  the  sun  if  these  relative  positions  and  the  corresponding  relative  velocities 
are  known  at  any  one  time.  Therefore,  what  we  ideally  want  in  GRT  is  differential  equations 
which  give  us  the  acceleration  of  each  body  (specified  by  some  parameters)  in  terms  of  the  rela¬ 
tive  positions  and  velocities  of  the  other  bodies.  By  the  very  nature  of  a  relativistic  theory  such 
hopes  are  doomed  except  in  the  limit  in  which  the  velocity  of  light  is  infinite,  i.e.,  in  the  "lowest” 
approximation.  There  are  three  reasons  for  this: 

(a)  Retardation.  Interactions  between  particles  take  place  via  fields  which 
propagate  at  the  speed  of  light.  Thus  even  if  forces  can  be  written 

to  depend  on  interparticle  distances  and  velocities,  they  must  refer  to 
retarded  quantities,  i.e.,  to  the  values  occurring  at  an  earlier  time 
such  that  signals  emitted  then  will  arrive  at  the  present  moment  at  the 
particle  of  interest.  In  order  to  express  everything  in  terms  of  vari¬ 
ables  at  a  single  time,  it  is  necessary  to  make  a  power  series  expansion 
in  the  retardation  time  and  cut  it  off  before  higher  derivatives  than  the 
velocity  appear.  There  are  further  complications  since  at  any  given 
time  some  of  the  energy  and  momentum  of  the  system  is  "in  transit"  so 
that,  for  example,  the  momentum  of  the  system,  which  is  approximately 
conserved,  is  not  just  the  sum  of  the  momenta  of  the  particles,  and  the 
expression  for  the  center  of  energy  (the  analogue  of  the  center  of  mass) 
becomes  much  more  complicated  than  is  the  case  in  Newtonian  theory. 

(b)  Proper  Time.  Relativistic  equations  of  motion  such  as  the  geodesic 
equation  (2-14),  always  involve  derivatives  with  respect  to  proper  time 
along  the  world  line  of  the  particle.  To  obtain  differential  equations  in 
terms  of  coordinate  time,  ds/dt  must  be  expanded.  This  involves  fix¬ 
ing  a  coordinate  system  and  the  resulting  expressions  will  involve  the 
"absolute"  velocities  of  all  the  particles  in  this  system  rather  than  just 
the  relative  velocities  of  the  particles.  There  seems  to  be  no  way  of 
avoiding  this.  In  the  standard  treatment  of  the  GRT  advance  of  perihe¬ 
lion,  the  fields  corresponding  to  the  sun  at  rest  are  used  so  the  absolute 
velocities  of  planets  become  velocities  relative  to  the  sun.  This,  however, 
ignores  the  forces  exerted  by  the  planets.  If  we  try  to  include  these 
forces  then  the  metric  becomes  time-dependent  and  we  can  no  longer 
exactly  identify  coordinate  time  with  proper  time  at  the  sun,  i.e.,  the 
"time"  in  the  static  Schwarzschild  solution,  and  so  we  run  into  the  same 
difficulties. 

(c)  Radiation.  This  is  a  basically  unsolved  problem.  If  we  have  a  system 
of  charged  particles  then  we  expect  them  to  radiate  so  that  some  kind  of 
"radiation  reaction"  force  must  be  included  in  addition  to  the  interparticle 
forces.  Since  the  radiation  reaction  force  involves  the  time  derivative 

of  the  acceleration,  it  radically  changes  our  concept  of  equations  of  mo¬ 
tion  and  in  fact  converts  them  from  differential  to  integro-differential 
equations  whose  properties  are  not  too  well  understood.30  Fortunately, 
since  we  are  forced  for  retardation  reasons  to  expand  in  powers  of  1 /c 
anyhow,  we  can  avoid  the  radiation  reaction  force  in  the  first  two  approx¬ 
imations31  because  the  radiation  reaction  force  is  proportional  to  1/c3. 

Thus,  if  we  are  careful  to  keep  only  terms  up  to  1/c2  we  can  "forget" 
about  radiation.  In  electrodynamics,  in  the  lowest  approximation  we  ob¬ 
tain  the  Coulomb  interactions  only,  while  in  the  second  approximation  we 
get  the  so-called  Darwin  Hamiltonian,  run  into  absolute  velocities,  and 
in  fact  get  very  similar  results  to  those  we  will  obtain  in  the  post-Newtonian 
approximation  to  GRT. 

A  careful  discussion  of  all  this,  including  an  analysis  of  the  Lorentz 
covariance  of  the  resulting  theory  (it  is  covariant  with  respect  to 
Lorentz  transformations  if  we  formally  expand  the  transformation  equa¬ 
tions  in  powers  of  l/c  and  keep  terms  only  up  to  1/c2)  can  be  found  in 
Fock’s  monograph.3 2  Sweeping  radiation  under  the  rug  in  this  manner 
is  fine  formally,  but  a  skeptic  might  ask  why  one  believes  the  approxi¬ 
mate  theory  bears  any  resemblance  to  the  full  theory.  The  only  reply 
one  can  make  is  to  express  the  pious  hope  that  radiation  by  slowly  moving 
charges  is  a  small  effect  and  can  be  ignored.  Whether  or  not  the  Darwin 
Hamiltonian  is  a  good  approximation  to  the  full  theory,  its  quantum 
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theoretical  analogue,  the  Breit  Hamiltonian,  seems  to  work  very  well 
for  two-electron  systems .33  The  situation  in  GRT  is  even  worse  since 
no  one  knows  what  radiation  means  or  whether  or  not  it  makes  sense  to 
talk  about  radiation  of  energy  and  momentum  by  an  isolated  system .22 
If  one  sticks  to  the  linearized  theory  there  is  a  complete  analogy  to 
electrodynamics  and  radiation  is  well  defined  and  goes  as  l/c^  (Ref.  34) 
so  one  can  again  sweep  it  under  the  rug  in  the  lowest  few  approxima¬ 
tions.  The  difficulty  is  that  no  one  has  any  idea  of  the  amount  of  corre¬ 
spondence  between  the  linearized  and  full  theories.  In  any  case,  we 
can  formally  obtain  the  post-Newtonian  equations  of  motion  without 
running  into  trouble,  but  it  is  difficult  to  go  beyond.3 5 

We  now  discuss  the  various  derivations  of  the  post-Newtonian  equations  of  motion  that 
appear  in  the  literature.  All  end  up  with  the  same  equations  although  the  derivations  are  super¬ 
ficially  different.  The  common  idea  is  to  expand  the  field  equations  in  powers  of  l/c  as  we  did 
in  deriving  the  Newtonian  equations.  Then  it  is  assumed  that  the  energy-momentum  tensor 
(1)  looks  like  a  sum  of  6-functions  corresponding  to  each  particle  or  (2)  represents  incoherent 
matter  with  the  matter  "concentrated”  near  the  location  of  each  particle  or  (3)  vanishes  outside 
of  a  number  of  nonintersecting  arbitrarily  narrow  tubes  in  space-time  within  which  it  becomes 
singular.  Much  verbal  blood  has  been  shed  by  proponents  of  the  different  energy-momentum 
tensors  but  it  is  obvious  that  they  are  all  rather  equivalent  and  so  it  is  not  surprising  that  they 
all  get  the  same  equations.  W'hether  one  uses  6-functions  or  "concentrated"  matter  is  obviously 
irrelevant  except  that  with  6-functions  one  has  to  drop  some  formally  infinite  terms  which  do 
not  appear  when  the  matter  is  merely  "concentrated."  The  derivation  based  on  solving  the  empty- 
space  equations  in  a  multiply-connected  region  outside  the  world-lines  of  the  sources  is  quite 
different  in  form. 

Consider  first  the  6-function  approach.20  We  have  seen  that  to  the  lowest  order  the  only 
component  of  the  metric  tensor  which  is  determined  by  comparison  with  Newtonian  theory  is 

g44=-l-^f  [Eq.  (2-53)] 

c 

3  6 

where  (p  is  the  Newtonian  potential.  The  other  components  are  highly  arbitrary'  since  if  we 
are  given 


gij 


T) . .  +  Xh. . 

l2 


we  can  perform  an  infinitesimal  coordinate  change 


[Eq.  (2-45)] 


— i  i  ,  , ,  i/  v 
x  =  x  +  X£  (x) 

with  £l(x)  constrained  only  to  vanish  at  infinity  and  to  satisfy 


(2-60) 
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which  is  necessary  to  preserve 
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(2-61) 


(2-62) 


Thus  both  the  barred  and  unbarred  coordinate  systems  are  equally  suitable  in  the  Newtonian 
limit,  but  the  metric  tensor  in  the  barred  system  has  components 
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g. .  =  g. .  —  X 

(fii  .  ffj) 

'  3x^  ax1 ' 

+  0(X2)  =  j,  +  X  hj. - 

(!ii  ♦  fill 

'axJ  ax1'! 

(2-63) 

=  T)..^ 

S1  'if 

(2-64) 

Because  of  the  restriction  of  Eq.  (2- 

61), 

g44  =  ^44  +  °(X  ] 


(2-65) 


but  all  the  other  h^  are  arbitrary  although  one  must  obviously  choose  some  values  for  them  in 

order  to  pin  down  the  coordinate  system  to  first  order.  It  is  conventional  to  choose  h..  to  make 

20 

the  first  approximation  to  the  metric 


ds2  =  (l  +  c2dt2  -  (l  -  j(dx2  +  dy2  +  dz" 


(2-66) 


The  historical  basis  for  this  choice  of  metric  is  that  the  equations  of  the  linearized  theory 
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of  gravitation  then  assume  their  simplest  form.  This  choice  is  almost  completely  analogous 
to  that  of  the  Lorentz  Gauge  for  the  potentials  in  electrodynamics,  the  coordinates  in  Eq.(2-66) 
being  known  as  "harmonic"  or  "de  Donder"  coordinates.  Except  for  Fock's  school,1^'17  few 
attach  any  special  significance  to  harmonic  coordinates  in  the  full  theory  (Fock's  point  of 
view  is  discussed  on  page  15) 
the  second  can  be  developed. 


view  is  discussed  on  page  15).  In  any  case,  once  the  first  approximation  to  h. .  has  been  chosen, 


Accordingly,  we  substitute  the  g.^  of  Eq.(2-66)  into 
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[Eq.  (2-42)1 


which  yields,  for  example,^® 


=  Jj  mac"”  l  1  + 


<5(r  —  R  ) 
a 


(2-67) 


and  solve  the  resulting  field  equations.  One  solution  correct  to  second  order  is 


.  20 


Zap  =  + 


(2-68) 
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(2-70) 


In  order  to  obtain  the  last  three  equations  it  is  necessary  to  ignore  all  terms  which  become 
infinite  at  the  positions  of  the  particles.  This  can  be  "justified"  by  looking  at  "concentrated" 
matter  instead  of  at  "point  masses."  The  equations  of  motion  of  a  test  particle  which  moves  in 
such  a  field  can  now  be  obtained;  they  are  just  the  equations  of  a  geodesic  and  therefore  follow 
from  the  variational  principle 
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<5  \  ds  =  0 


[  Eq.  (2-11)] 


which  we  may  rewrite  to  introduce  a  Lagrangian 
6  ^  Ldt  =  0 


(2-71) 


T  ds 

L="mc  dT 


(2-72) 


where  the  coefficient  "-me"  is  conventional.  Although  it  is  not  necessary  for  our  purposes,  one 
can  obtain  a  Lagrangian  for  the  N-body  system.  This  is  not  the  sum  of  the  Lagrangians  for  each 
body  but  is  constructed  so  that  it  leads  to  the  correct  force  on  each  body  for  a  given  motion  of 
the  other  bodies.  The  result  is 
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(2-74) 


and  all  terms  in  the  sum  which  blow  up  are  understood  to  be  omitted. 

In  this  way  the  GRT  N-body  problem  can  be  reduced  in  second  approximation  to  a  problem 

in  classical  mechanics  with  the  Lagrangian  in  Eq.(2-73)  replacing  the  one  of  ordinary  Newtonian 

theory.  The  Lagrangian  belongs  to  the  standard  class  of  those  which  involve  velocity-dependent 

forces  and  a  Hamiltonian  can  be  constructed,  conservation  laws  derived,  etc.,  in  the  usual  way. 

Since  we  are  just  interested  in  the  equations  of  motion,  we  will  not  discuss  these  points  here, 
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but  an  exhaustive  analysis  may  be  found  in  the  monograph  of  Infeld  and  Plebanski. 

There  is,  however,  one  very  important  result  which  follows  from  the  fact  that  the  post- 
Newtonian  equations  of  motion  are  derivable  from  the  Lagrangian  Eq.(2-73).  In  the  discussion 
of  the  Newtonian  approximation  we  saw  that  although  the  equations  of  motion  were  unique,  the 
coordinate  system  was  not.  A  similar  result  is  true  in  the  post-Newtonian  approximation.  The 
post-Newtonian  g.  .  [Eqs.  (2-68)  to  (2-70)]  are  far  from  unique  even  with  the  asymptotic  boundary 
condition;  the  general  solution  involves  four  arbitrary  functions  which  go  to  zero  sufficiently 
rapidly  at  spatial  infinity.  If  one  computes  the  Lagrangian  which  replaces  Eq.  (2-63)  when  the 

general  solution  is  used,  however,  one  finds  that  it  differs  from  Eq.(2-73)  by  a  total  time 
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derivative,  and  it  is  well  known  from  classical  mechanics  that  this  means  that  both  Lagrangians 
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lead  to  the  same  equations  of  motion.  It  is  not  known  whether  similar  results  hold  in  higher 
orders  because  of  the  ambiguities  involving  gravitational  radiation. 

Before  discussing  the  equations  of  motion  themselves,  let  us  examine  in  some  detail  the 
other  derivations  of  the  equations  of  motion  which  appear  in  the  literature.  We  will  first  outline 
Fock's  derivation,1  ^  which  is  straightforward,  and  then  his  motivations,  which  are  controversial.40 
Fock  starts  with  the  field  equations 
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(Eq.  (2-1)1 


D  _  —  a  R  -  T 

ik  2  ^ik  ~  4  ^ik 

c 

which  by  the  definition  of  R  in  Eq.  (2-24)  can  be  written  as 


R  _  87tG 

Rik  ’  ~ 
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,Tik  -  I  elkT>  • 


(2-75) 


He  now  decomposes  R^  into  a  sum  of  terms,  none  of  which  separately  transform  as  tensors. 
Thus,  symbolically. 


Eik  *  Aik  +  Bik  =  <T,k  -  5  *ikTI  • 


(2-76) 


Fock  finds  that  if  a  particular  class  of  coordinate  systems,  the  so-called  "harmonic”  or 

"de  Donder,"  systems,  are  used,  then  the  in  Eq.(2-76)  vanish.  Thus  the  Einstein  equations 

can  be  replaced  by  the  equations 


Aik  "  4  (Tik  2  gikT) 
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C.  =  0 
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(2-77) 


(2-78) 


where  Eqs.  (2-78)  are  four  equations  which  must  be  satisfied  by  g. .  in  order  to  make  the  metric 
"harmonic."  Equations  (2-77)  and  (2-78)  are  no  longer  "tensor"  equations  since  A.^  is  not  a 
tensor,  but  are  still  a  perfectly  respectable  set  of  partial  differential  equations.  There  is  no 
longer  a  need  to  worry  about  matters  such  as  the  Bianchi  identity  when  solving  Eq.  (2-77)  but  of 
course  a  solution  of  Eq.  (2-77)  is  no  longer  necessarily  a  solution  of  the  Einstein  field  equations; 
it  is  only  a  solution  if  it  also  satisfies  Eq.  (2-78),  in  which  case  of  course  it  will  satisfy  the 
Bianchi  identity  in  the  harmonic  coordinate  system  since  A.^  will  be  the  Ricci  tensor.  As  usual, 
all  this  has  an  analogue  in  classical  electrodynamics.  Maxwell's  equations  can  be  written  in 
terms  of  the  four  potentials.  Given  any  set  of  potentials  which  satisfy  the  equations,  any  other 
potentials  related  to  the  originals  by  a  gauge  transformation  yield  the  same  fields.  If  an  auxil¬ 
iary  condition  is  imposed  on  the  potentials  (e.g.,  the  Lorentz  condition)  so  as  to  simplify  the 
equations,  then  the  resulting  equations  for  the  potentials  are,  in  the  first  place,  not  completely 
gauge  invariant.  In  the  second  place,  a  solution  of  these  equations  is  not  in  general  a  solution 
of  Maxwell's  equations  but  only  becomes  so  if  it  satisfies  the  auxiliary  condition  as  well. 

Roughly  speaking,  what  Fock  has  done  is  to  choose  his  coordinate  condition  to  all  orders, 
ab  initio,  rather  than  to  choose  it  in  each  order  as  has  been  done  in  the  derivations  of  the  New¬ 
tonian  and  post-Newtonian  equations  outlined  above.  By  solving  Fock's  equations  in  the  Newtonian 
and  post-Newtonian  approximations,  the  end  result  is,  as  mentioned  before,  the  exact  same 
equations  of  motion. 

Fock  believes  that  Einstein's  field  equations  are  incomplete  and  that  they  should  be  supple¬ 
mented  with  the  requirement  that  only  a  special  class  of  harmonic  (or,  in  an  iterative  solution, 
approximately  harmonic)  coordinates  be  used  because  they  are  the  most  "physical."  In  order  to 
clarify  this,  consider  the  case  of  a  flat  space-time.  Suppose  we  decide  to  use  spherical  coordi¬ 
nates  in  place  of  Cartesian  coordinates.  Then  the  mathematics  singles  out  the  point  at  the  origin 
of  the  spherical  coordinate  system  (e.g.,  the  transformation  from  Cartesian  to  spherical  coordi¬ 
nates  is  singular  at  the  origin  since  0  and  cp  are  not  well  defined  there)  even  though  it  is 
"physically"  just  another  point  in  space  time  and  geometrically  indistinguishable  from  any  other 
point. 
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Mathematically  speaking,  a  flat  space-time  has  a  high  degree  of  symmetry  since  every 


in  Cartesian  coordinates  reflects  this  more  clearly  while  any  other  coordinates  will  introduce 
’’unphysical"  singularities  into  the  metric.  This  does  not  usually  hamper  physicists  calculating 
in  spherical  coordinates  since  the  singular  properties  of  the  origin  can  more  or  less  automat¬ 


ically  be  taken  into  account  by  modifying  boundary  conditions  (e.g.,  in  quantum  mechanics  the 


wave  function  in  Cartesian  coordinates  must  be  finite  everywhere  whereas  solutions  of  the  Dirac 
equation  in  spherical  coordinates  are  allowed  to  be  mildly  singular  at  the  origin).  It  is,  however, 
more  serious  in  GRT  where  the  metric  itself  is  being  studied  to  separate  "physical”  from 
"coordinate"  singularities. 

It  is  found  that  if  Fock's  equations  are  solved  for  an  empty  space-time  which  is  Minkowskian 
at  spatial  infinity,  the  resulting  coordinates  are  just  a  set  of  Cartesian  coordinates  (or,  to  be 
more  precise,  the  family  of  Cartesian  coordinates  related  to  each  other  by  inhomogeneous 
Lorentz  transformations).  A  similar  theorem  can  be  proved  even  if  there  is  matter  present 
providing  a  static  solution  (i.e.,  one  for  which  the  g. .  have  no  explicit  time  dependence)  exists 
in  which  case  Fock's  coordinates  become  asymptotically  Cartesian  at  spatial  infinity.  Fock 
wishes  to  generalize  this  and  in  particular  to  use  only  harmonic  coordinates  in  discussing  the 
problem  of  gravitational  radiation.  It  is  at  this  point  that  many  other  relativists  part  company 
with  Fock,  for  if  anything  like  gravitational  radiation  exists  then  there  is  no  compelling  reason 
to  assume  that  space-time  at  infinity  is  asymptotically  flat  and  therefore  there  is  no  compelling 
reason  to  consider  solutions  in  harmonic  coordinates  as  being  freer  of  spurious  "coordinate 
singularities"  than  others.  Accordingly,  most  relativists  feel  that  although  harmonic  coordinates 
may  be  convenient  in  particular  problems  where  one  is  looking  for  a  solution  which  is  of  the  form 
of  a  static  solution  plus  a  small  perturbation,  they  do  not  have  any  special  significance  in  more 
general  cases. 
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Finally,  there  is  the  Einstein-Infeld-Hoffman  (EIH)  method.  ’  '  Like  the  other  methods 

described  above,  it  is  based  on  a  power  series  expansion  in  l/c,  the  assumption  that  the  metric 
is  asymptotically  Minkowskian  at  spatial  infinity,  and  the  assumption  that  all  time  dependence 
is  due  to  the  motion  of  the  particles.  In  its  original  formulation,  the  world-lines  of  the  particles 
were  assumed  not  to  intersect,  were  each  enclosed  in  a  "tube"  and  the  field  equations  over  all 
space-time  were  then  replaced  by  the  equations,  with  no  matter  present,  extended  over  the 
multiply-connected  region  of  space-time  outside  the  tubes. 

One  solution  is  obviously  a  flat  space-time,  i.e.,  one  in  which  nothing  is  singular  anywhere. 
Since  this  is  what  we  expect  if  there  is  no  matter  at  all,  this  solution  is  rejected  as  unphysical 
and  solutions  are  sought  which  become  singular  on  the  world-lines  of  the  particles.  The  sur¬ 
prising  feature  of  GRT  is  that  there  are  no  solutions  of  the  equations  for  arbitrary  singularities 

and  world-lines.  This  is  a  consequence  of  the  Bianchi  identities.  To  see  how  this  comes  about, 
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we  again  turn  to  an  analogy  from  electrodynamics.  Suppose  that  the  free-field  Maxwell  equa¬ 
tions  hold  everywhere  outside  the  world-lines  of  charged  point  particles  and  that  these  world-lines 
do  not  intersect.  The  charge  on  one  of  these  particles  at  a  certain  time  can  then  be  defined  by 
Gauss'  law  as 


(2-79) 
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where  the  integral  is  extended  over  a  closed  spatial  surface  surrounding  the  particle  at  this 
time  and  does  not  surround  any  other  particles  (this  can  be  done  because  of  the  nonintersection 
hypothesis).  Then,  making  use  of  Maxwell’s  equations  in  free  space  and  Stokes’  theorem  yields 


dq  _  1  fa  E 
dt  "  4tt  ?  at 


V  x  B  •  dS  =  0 


(2-80) 


Thus  there  can  be  no  solutions  to  Maxwell’s  equations  in  which  the  charges  of  isolated  particles 
change.  This  is  of  course  just  a  special  case  of  the  well-known  result  that  one  can  derive  the 
equation  of  continuity  from  Maxwell’s  equations  and  that  therefore  specified  charge  and  current 
distributions  must  be  consistent  with  the  equation  of  continuity  or  else  there  are  no  solutions  to 
Maxwell's  equations. 

There  is  one  further  important  point  to  make.  Suppose  the  particle  does  not  carry  a  net 
charge  but  is  a  dipole.  We  can  visualize  this  dipole  as  consisting  of  two  equal  and  opposite 
charges  which  are  so  ’’close  to  each  other"  that  the  surface  surrounding  the  particle  always  in¬ 
cludes  both  charges.  Then  the  integrals  in  Eqs.  (2-79)  and  (2-80)  are  always  zero  and  solutions 
can  be  found  to  Maxwell’s  equations  in  which  the  strength  of  the  singularity,  i.e.,  its  dipole 
moment,  varies  in  time  (imagine  the  two  charges  exchanging  charge).  So  the  equation  of  con¬ 
tinuity  by  itself  imposes  a  restriction  only  on  a  part  of  the  possible  singularity.  Since  the 
electric  field  of  a  point  charge  is  spherically  symmetric  and  is  an  inverse  square  law  force 
while  the  field  of  a  dipole  or  higher  multipole  is  not  spherically  symmetric  and  is  an  inverse 
cube  (or  higher)  law  force,  the  nature  of  the  singularity  can  be  identified  by  looking  at  how  sin¬ 
gular  E  becomes  as  the  size  of  the  tube  decreases  to  zero.  The  multipole  structure  of  the  sin¬ 
gularity  can  be  "identified,"  and  it  can  be  said  that,  e.g.,  only  a  charge  or  only  a  dipole  exists. 

Obviously,  this  analysis  depends  critically  on  the  linearity  of  Maxwell’s  equations.  The 
analogues  of  the  equation  of  continuity  in  GRT  are  the  four  Bianchi  identities  which  act  as  con¬ 
sistency  conditions  and  certainly  impose  four  constraints  on  the  world-lines  of  singularities. 

The  fact  that  the  theory  is  nonlinear  and  has  the  same  form  in  any  coordinate  system  wreaks 
havoc  with  the  multipole  expansion;  in  fact,  the  identification  of  singularities  by  means  of  the 
dependence  of  the  fields  on  the  "size"  of  the  world  tube  as  it  shrinks  to  a  line  can  only  be  done 

in  low  orders  of  approximation  and  it  has  so  far  proved  impossible  to  characterize  the  properties 
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of  mass  distributions  in  an  invariant  manner.  A  not  entirely  independent  problem  that  also 
arises  is  directly  due  to  the  nonlinearity  of  the  equations  which  results  in  singularities  being 
piled  on  top  of  singularities  at  each  level  of  approximation.  These  can  be  argued  away  in  the 
lowest  orders  by  an  averaging  process  (called  "tweedling"  by  Infeld)13,45  but  again  this  is  a  non¬ 
invariant  procedure. 

In  the  original  EIH  paper,  no  assumptions  were  made  at  the  outset  about  the  nature  of  the 

singularities,  but  it  was  found,  not  unexpectedly,  that  assumptions  had  to  be  made  about  how- 
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singular  the  fields  would  be  on  the  surface  of  the  tube  and  EIH  assumed  the  least  nontrivial 
1 1 

singularity.  Infeld  showed  that  this  was  equivalent,  at  least  up  to  the  post-Newtonial  approxi¬ 
mation,  to  assuming  a  6-function  distribution  of  matter  of  the  form  of  Eq.  (2-41)  and  that  by  using 
Eq.  ( 2-41 )  from  the  beginning  the  calculation  could  be  simplified  tremendously. 

We  will  now  outline  the  essential  features  of  Infeld' s  calculation.  An  expansion  is  made  of 
g. .  as  in  Eq.  (2-43)  with  l/c  as  the  parameter  of  smallness,  and  of  the  energy  momentum  tensor 
of  Eq.  (2-41).  If  all  the  terms  of  a  given  order  are  then  collected,  the  field  equations  are  re¬ 
placed  by  an  infinite  set  of  coupled  equations  each  of  which  has  the  semi-symbolic  form 
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(2-81) 
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+  N(gim)  +  aTim  =  0 


Both  the  quantity  N(gim),  which  represents  the  terms  of  a  given  order  which  arise  from  multi¬ 
plying  terms  of  lower  order  as  required  by  the  nonlinearity  of  the  equations,  and  the  terms  T 
which  come  from  the  energy  momentum  tensor,  involve  only  quantities  of  lower  order.  Thus 
the  only  new  h. .  appear  in  the  first  term.  K  is  linear  in  these  quantities  and  in  fact  looks  just 
like  the  linear  part  of  the  Riemann  tensor  of  Eq.  (2-22).  This  is  very  useful  for  knowing  just 
how  much  freedom  there  is  in  choosing  a  coordinate  system,  up  to  this  order,  which  will  simplify 
the  equations.  It  also  follows  that  Eq.  (2-81)  is  consistent  only  if  N  +  aT  satisfies  Bianchi-like 
identities  with  K  regarded  as  a  Riemann-like  tensor.  This  leads  to  constraints  on  the  motion  of 
the  sources  which  turn  out  to  be,  successively,  the  Newtonian  and  post- Newtonian  equations.  In 
higher  orders,  difficulties  arise  because  some  of  the  terms  in  N  do  not  go  to  zero  fast  enough 
at  spatial  infinity  and  so  the  whole  program  must  be  re-examined. 

To  summarize  the  state  of  the  art,  the  physically  reasonable  assumptions  that  the  gravita¬ 
tional  field  produced  by  a  system  of  structureless  point  particles  should  vanish  at  infinity,  have 
a  time  dependence  only  through  the  time  dependence  of  the  positions  of  the  particles,  and  become 
singular  at  the  position  of  the  particle  in  a  "spherically  symmetrical"  manner,  lead  to  unambig¬ 
uous  equations  of  motion  up  to  and  including  the  post-Newtonian  approximation  (if  one  makes  the 
additional  assumption  that  the  kinetic  energy  of  each  particle  is  of  the  same  order  of  magnitude 
as  the  Newtonian  potential  energy,  this  last  assumption  leads  to  equations  of  motion  which  are 
symmetric  in  each  of  the  particles).  It  is  not  at  all  clear  to  what  extent  this  procedure  can  be 
extended  beyond  the  post-Newtonian  approximation  or  whether  it  is  in  any  sense  a  "good"  approxi¬ 
mation  to  the  exact  solution  from  a  mathematical  point  of  view.  There  are  further  difficulties 

.  .  44 

involved  in  extending  the  analysis  to  include  particles  with  structure,  e.g.,  spinning  particles. 
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ni.  THE  POST- NEWTONIAN  EQUATIONS  OF  MOTION 


The  post-Newtonian  equations  of  motion  can  be  derived  in  the  standard  manner  from  the 
Lagrangian  of  Eq.  (2-73)  by  calculating  the  left-hand  side  of  the  equation 


where  X  and  V  are  respectively  the  position  and  velocity  vectors  of  the  a^  particle  and  where 
a  a 


d  def  9f 

®f<x .-v  ■  z* 

a  a 


V  + 
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8f 

dV 


dV 

_ c 

dt 


(3-2) 


The  Eqs.  (3-1)  will  differ  markedly  in  appearance  from  the  Newtonian  equations.  Since  the 
Lagrangian  contains  no  higher  order  time  derivatives  than  the  first,  the  equations  will  still  con¬ 
tain  no  higher  order  derivatives  than  the  second  and  will,  in  fact,  be  linear  in  the  accelerations. 
However,  because  of  the  complicated  dependence  of  the  Lagrangian  on  the  positions  and  velocities 
of  the  particles  and  the  fact  that  the  time  derivative  is  to  be  calculated  according  to  Eq.  (3-2), 
each  equation  of  the  set  will  involve  the  acceleration  of  all  of  the  particles.  After  a  somewhat 
tedious  calculation,  the  following  equations  of  motion  are  obtained: 
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and  the  sums  are  extended  over  all  the  particles  but  the  ath.  We  have  further  simplified  the 
equations  by  using  units  in  which  the  speed  of  light  and  the  gravitational  constant  are  unity, 
i.e.,  in  Eq.  (3-3), 

"v"  =  v/c  (3-5) 


and 


=  GM/c2 


(3-6) 


We  are  not  really  interested  in  the  general  solution  of  Eq.  (3-3)  but  only  in  one  that  looks 

like  a  perturbed  Newtonian  solution.  We  are  therefore  justified  in  substituting  the  Newtonian 

values  of  the  acceleration  in  those  terms  on  the  right-hand  side  of  Eq.  (3-3)  which  contain  the 

accelerations  explicitly,  since  the  accelerations  are  multiplied  by  factors  of  the  formal  order 
2  2 

of  v  /c  or  the  Newtonian  gravitational  potential,  and,  from  the  derivation  of  the  post-Newtonian 
approximation  as  discussed  in  Sec.  I,  any  errors  incurred  by  this  substitution  are  then  of  the 
same  formal  order  as  terms  which  have  already  been  neglected.  If  we  substitute 
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on  the  right-hand  side  and  reorder  the  terms  slightly,  we  obtain 
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(3-8) 


Substantially  the  same 


Equations  (3-8)  were  first  obtained  by  Eddington  and  Clark  in  1938. 

47 

equations  had  been  obtained  in  1915  by  Droste  and  De  Sitter  who  unfortunately  made  a  mistake 
and  so  obtained  the  wrong  coefficients  for  the  last  two  terms  in  Eqs.  (3-8).  Wc  will  apply  these 
equations  to  find  the  GRT  corrections  to  the  planetary  equations  of  motion  and  to  the  earth-moon- 
sun  system. 
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IV.  GRT  AND  THE  MOTION  OF  THE  PLANETS 


Although  the  equations  of  motion  are  symmetric  in  each  of  the  bodies,  there  is  an  inherent 
asymmetry  in  the  solar  system  due  to  the  fact  that  the  sun  is  so  much  more  massive  than  the 
combined  mass  of  all  the  other  bodies.  Thus,  terms  of  the  same  formal  order  in  l/c  can  have 
vastly  different  numerical  values  depending  on  whether  a  solar  or  planetary  mass  is  involved. 
As  a  first  step,  let  us  rewrite  Eqs.  (3-8)  for  the  i^  planet  to  bring  out  the  dependence  on  the 
solar  mass. 
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where  now  the  subscript  "oM  refers  to  the  sun  and  the  solar  term  is  excluded  from  all  sums. 
A  similar  equation  can  be  written  for  the  sun: 
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where  the  sums  run  over  all  the  planets. 

Now,  Eqs.  (4-1)  and  (4-2)  involve  more  quantities  than  we  can  measure  since  we  must  choose 
an  origin  for  our  coordinate  system.  If  we  look  at  the  Newtonian  terms  alone,  i.e.,  ignore  the 
right-hand  side  of  Eqs.  (4-1)  and  (4-2),  and  subtract  Eq.  (4-2)  from  (4-1)  we  obtain  the  system  of 
equations 
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dt 
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Thus,  in  the  Newtonian  approximation,  we.  can  obtain  a  set  of  equations,  the  ’’planetary  equations," 
involving  only  relative  quantities.  Furthermore,  by  multiplying  each  Eq.  (4-1)  by  m.  and  Eq.  (4-2) 
by  M  and  adding  all  the  equations  together  in  the  Newtonian  approximation,  we  obtain 


y  m.a.  +  Ma  =  Y  m.(iT.  —  a)+(M+ym.  )a  =0 

L  1  1  °  u  1  1  °  \  “  7  ° 


(4-4) 


or 


The  last  equation  states  that  the  Newtonian  center  of  mass  is  unaccelerated.  Eqs.  (4-3)  and  (4-5) 
are  equivalent  to  Eqs.  (4-1)  and  (4-2)  in  the  Newtonian  approximation.  The  center  of  mass  equa¬ 
tion,  unlike  the  planetary  equations  (4-3),  does  not  involve  only  relative  quantities. 

The  planetary  equations  do  not  involve  the  velocity  of  the  center  of  mass  and  look  like  a  com¬ 
plete  self-contained  set  of  equations.  This  is  illusory.  The  masses  which  appear  are  parameters 
which  must  be  determined  by  fitting  the  observed  motion  to  the  predicted  motion  with  arbitrary 
masses;  to  do  this  one  must  have  a  theory  of  light  and  of  other  electromagnetic  phenomena. 
Newtonian  mechanics  does  not  include  anything  about  light  so  additional  postulates  are  needed. 

The  usual  postulate  is  that  light  travels  in  a  straight  line  with  the  same  speed  ”c"  as  is  obtained 
in  terrestrial  measurements  of  the  speed  of  light.  Once  this  assumption  is  made,  the  velocity  of 
the  center  of  mass  indeed  becomes  irrelevant. 
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When  we  look  at  the  GRT  equations,  we  see  that  although  only  relative  positions  appear, 
absolute  velocities  are  involved  (this  is  inherent  in  higher  non-relativistic  approximations  to  a 
relativistic  theory  as  was  discussed  in  Sec.  I).  Now,  we  know  that  ignoring  the  relativistic  terms 
completely  allows  one  to  obtain  an  excellent  fit  to  the  motion  of  the  planets  around  the  sun  and 
that  the  principal  remaining  discrepancies  can  be  removed  by  adding  on  the  GRT  Schwarzschild 
field  produced  by  the  sun,  treating  the  sun  as  if  it  were  at  rest.  There  will  thus  be  a  very  good 
solution  in  which  the  "absolute”  velocity  of  the  sun  is  much  smaller  than  the  velocities  of  the 
planets  relative  to  the  sun.  We  will  make  the  assumption  that  if 

def 

w  =  "absolute  velocity"  of  the  sun  (4-6) 

then 

(Mgun)  (w)  <  (Mpjanet)  ^ve^oc^y  planet  relative  to  sun)  .  (4-7) 

This  assumption  cannot  be  rigorously  justified:  an  attempt  can  only  be  made  to  try  to  see 
whether  the  equations  derived  using  the  assumption  describe  the  observed  motions  of  the  planets. 
It  is  known  that  w  can  be  set  equal  to  zero  and  good  results  still  obtained,  and  that  the  estimate 
we  use  is  compatible  with  this.  Another  heuristic  reason  for  this  estimate  is  that  it  can  give 
the  Newtonian  center  of  mass  a  small  velocity  as  can  be  seen  from  Eq.  (4-5).  Equation  (4-5)  is 
only  approximately  true  when  the  full  post-Newtonian  equations  are  used,  but  Eddington  and 
Clark  have  shown  that  there  is  no  secular  acceleration  of  the  Newtonian  center  of  mass  for  the 
two-body  problem  in  the  post- Newtonian  approximation  and  it  may  well  be  conjectured  that  this 
is  true  for  the  many-body  problem. 

In  the  case  of  several  bodies,  it  is  known  that  there  is  a  "center  of  energy"  which  moves 

with  constant  velocity  and  that  the  "center  of  energy"  reduces  to  the  Newtonian  center  of  mass 

when  we  go  from  the  post-Newtonian  to  the  Newtonian  approximation;48  Fock4^  shows  that  for  any 

reasonable  estimates  of  the  sizes  of  the  terms,  the  distinction  between  the  "center  of  energy" 

and  the  Newtonian  center  of  mass  can  be  ignored  in  the  two-body  problem.  We  will  assume  the 
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same  is  true  in  general.  This  is  also  the  traditional  assumption  made  by  De  Sitter  who  was 
the  first  to  look  at  the  implications  of  the  many-body  equations  for  the  solar  system.  De  Sitter 
even  made  the  further  assumption  that  the  "absolute  velocity"  of  the  Newtonian  center  of  mass 
vanished  since  he  wanted  the  most  "symmetric"  metric  at  large  distances  from  the  solar  system. 

If  the  above  discussion  of  "absolute  velocities"  and  the  like  seems  confused  and  not  very  con¬ 
vincing,  it  only  reflects  in  the  author's  opinion,  the  confused  and  not  very  convincing  state  of  the 
art.  To  obtain  our  equations  we  have  had  to  assume  that,  in  some  coordinate  system,  the  abso¬ 
lute  velocities  of  our  particles  are  small  compared  to  the  speed  of  light  and  we  have  set  things  up 
so  that  far  away  from  the  solar  system  space-time  becomes  flat.  But  of  course  if  we  go  far 
enough  away  we  come  to  other  stars  and  regions  where  space-time  is  no  longer  flat.  Thus  our 
coordinate  system  is  only  valid  in  a  region  which  is  large  compared  to  the  solar  system  but  small 
compared  to  interstellar  distances,  and  our  arguments  about  asymptotic  behavior  must  be  rein¬ 
terpreted  in  the  light  of  this  limitation.  It  may  be,  for  example,  that  the  freedom  to  choose  a 
coordinate  condition  in  each  order  is  not  really  there  but  that  there  exist  preferred  coordinate 
systems,  e.g.,  if  there  is  a  "center  of  mass"  of  the  universe,  a  coordinate  system  in  which  this 
point  is  fixed  may  be  the  one  to  use,  and  therefore  we  should  not  be  too  cavalier  in  our  treatment 
of  absolute  velocities  of  bodies  in  the  solar  system. 
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These  problems  of  course  are  not  peculiar  to  relativistic  astronomy;  even  in  Newtonian 
mechanics  the  planetary  equations  should  be  modified  if  the  solar  system  is  accelerating  relative 
to  the  "fixed  stars."  We  can  summarize  the  situation  as  follows:  We  have  obtained  equations  of 
motion  by  assuming  the  solar  system  contains  all  the  matter  in  the  universe.  In  reality  it  does 
not  and  no  one  knows  how  to  modify  the  theory  to  include  this  fact  or  to  what  extent  the  equations 
of  motion  and  conservation  laws  of  the  theory  are  affected.  So  we  are  not  breaking  the  rules  of 
the  game  by  making  any  "reasonable"  assumption  about  the  "absolute  velocity"  of  the  solar  sys¬ 
tem  and  looking  at  the  equations  which  result. 

Since  we  assume  the  Newtonian  center  of  mass  is  unaccelerated,  we  will  work  with  Eq.  (4-1) 
for  the  planetary  accelerations  and  replace  Eq.  (4-2)  for  the  acceleration  of  the  sun  by 
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If  Eqs.  (4-1)  are  written  as 
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which  defines  F.,  then  by  subtracting  Eq.  (4-2)  from  Eq.  (4-9)  we  obtain 
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We  know  we  want  a  solution  which  represents  a  perturbed  Newtonian  one.  We  therefore  re¬ 
place  Eq.  (4-10)  by  a  series  of  approximate  equations  obtained  by  collecting  terms  on  the  right 
which  are  of  the  same  order.  We  determine  the  order  of  a  term  by  noting  that  a  typical  velocity 
of  a  planet  relative  to  the  sun  is  about  30  km/sec.  Since  the  velocity  of  light  is  about  3X10 
km/sec,  this  gives  v.q/c  ~  10”4.  This  is  also  about  the  ratio  of  the  mass  of  a  planet  to  the  mass 
of  the  sun.  Finally,  the  velocities  relative  to  the  sun  satisfy  the  condition 
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which  follows  from  the  virial  theorem  for  an  inverse  square  force.  We  can  accordingly  classify 
accelerations  as  follows: 
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M 


mass  of  the  sun 


def 


m 

r 

v 

w 


def 

def 

def 

def 


mass  of  a  planet 

typical  interplanetary  or  planet-sun  distance 
typical  velocity  of  a  planet  relative  to  the  sun 
hypothetical  "absolute  velocity  of  the  sun"  [cf.  Eq.  (4-7)] 


(4-16) 


We  will  not  bother  going  beyond  the  third  order  and  the  author  even  has  some  doubts  about 
how  seriously  one  should  take  the  third  order  equations  since  w  may  appear  in  them.  Of  course 
the  w  terms  appear  in  third  order  only  because  we  have  guessed  the  size  of  w  to  be  given  by 
Eq.  (4-7),  and  if  this  is  wrong  they  could  first  appear  in  a  different  order.  We  will  see  that  the 
first  order  equations  are  the  usual  planetary  equations  [Eqs.(4-3)],  while  the  second  order 
corrections  give  the  terms  which  produce  the  relativistic  advance  of  the  perihelion  of  Mercury 
and  the  other  planets.  It  follows  that  the  influence  of  third  order  terms  will  probably  be  very 
hard  to  detect  and  fourth  order  terms  will  undoubtedly  be  irrelevant.  If  we  wish  to  include  third 
order  terms  we  would  have  to  assume  a  definite  value  for  w  and  the  only  assumption  which 
immediately  suggests  itself  is  De  Sitter’s,  viz.,  that  the  absolute  velocity  of  the  center  of  mass 
of  the  solar  system  vanishes.  This  determines  w,  which  then  has  the  magnitude  assumed  in 
Eq.  (4-7)  and  makes  the  third  order  terms  well-determined.  We  shall  write  the  third  order  equa 
tions  including  the  explicit  dependence  on  w. 


In  obtaining  the  equations  of  our  hierarchy  it  is  helpful  to  note  that  if  we  define 
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It  is  now  a  straightforward  task  to  identify  the  value  of  F.  up  to  a  given  order: 
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(4-23) 


Using  Eqs.  (4-18)  and  (4-19)  we  find,  in  successive  orders, 
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and 
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dt  _3  dt  _»2  \~3  —2'  all  j  X-2  -l' 


(4-27) 


These  are  our  final  equations.  Except  for  an  error  in  the  coefficients  in  the  third  order 
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acceleration,  they  were  first  found  by  Droste  and  De  Sitter. 
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V.  GRT  AND  THE  EARTH-MOON  SYSTEM 


Since  the  earth  and  moon  are  so  close,  we  cannot  use  the  order  of  magnitude  decomposition 
of  the  last  section  for  terms  involving  the  earth-moon  distance.  On  the  other  hand,  we  can  cer¬ 
tainly  ignore  third  order  terms  arising  from  the  sun  and  other  planets  in  the  equations  of  motion 
of  the  moon  relative  to  the  earth,  and,  if  we  desire,  include  them  in  the  equations  of  motion  of 
the  earth-moon  barycenter  by  merely  using  the  sum  of  the  masses  of  the  earth  and  moon  as  the 
"mass"  of  the  "barycenter."  Since  the  "absolute  velocity"  of  the  sun  does  not  appear,  presum¬ 
ably,  till  the  third  order,  we  may  as  well  choose  our  coordinate  system  centered  in  the  sun. 

We  define 
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X]VI  =  coorcRnates  °f  moon  relative  to  sun; 
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where  a 


(NR) 


»  /  p  \ 

is  the  nonrelativistic  part  and  a  the  GRT  part  of  the  acceleration  and  only  the 


-*■(  R ) 

sun,  moon,  and  earth  are  considered  in  a'  \ 

( R )  ( R) 

The  post-Newtonian  equations  give  us  expressions  for  a  ^  and  a^  .  Rather  than  writing 
these  down,  we  write  down  the  equivalent  expressions  for  the  relative  earth-moon  motion  and 
the  motion  of  the  earth-moon  barycenter: 
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(5-6) 


(3m  4-  4p) 
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For  the  earth-moon  barycenter  equations,  let 
del* 

mr^  =  m  +  }jl 


(5-7) 
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Essentially  these  equations  were  again  first  obtained  and  studied  by  De  Sitter  who  found 
a  very  interesting  effect,  essentially  due  to  the  velocity-dependence  of  the  Schwarzschild  force 
of  the  sun,  which  secularly  adds  an  extra  angular  velocity  to  earth  satellites,  e.g.,  the  moon 
We  will  discuss  this  effect  here  because,  on  the  basis  of  conversations  with  several  relativists, 
the  author  feels  that  knowledge  of  this  effect  is  not  widespread.  (The  only  text  which  discusses 
it  is  Eddington's.  He  calls  it  the  De  Sitter-Schouten  Effect.50)  To  see  the  effect  in  its  purest 
form,  let  us  look  at  the  terms  in  Eq.  (4-6)  which  are  due  to  the  sun  alone,  i  e.,  which  are  present 
if  we  set  the  masses  of  the  earth  and  of  the  moon  equal  to  zero.  Assume  also  that  the  earth 
moves  in  a  circle  about  the  sun  whose  radius  is  determined  by  Newton's  Law  for  circular  orbits: 


XE 

=  R 

(5-9) 

XM 

=  R  +  f 

(5-10) 

=  co  X  R 

(5-11) 
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-  M 

(5-12) 
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Idealizing  further,  let  us  look  for  relativistic  terms  which  survive  when  £  —  0f  i.e.,  we  are 
looking  for  GRT  effects  on  the  motion  of  a  very  close  satellite  orbiting  around  a  planet  moving 
in  a  circle. 

From  Eq.  (5-7)  and  Eqs.  (5-9)  through  (5-11)  we  then  find  that  the  only  terms  that  survive 

are 


d2r(R)  _  2m 


2(w  X  R) (R  •  +  R(R  •  u  x  $£) 


(5-13) 


dt  R 

Since  we  are  looking  for  secular  terms,  we  average  this  over  the  earth's  motion.  If  the 


plane  of  the  earth's  orbit  is  called  the  x-y  plane  then  obviously 

2 

<r  2>  =  <Ry2>  =  |«R  2>  +  <Ry2»  =  {<r2>  =  }Jr 

<RxRy>  =  <RxRz>  =  <RyRz>  =  <Rz>  =  0 
or,  in  matrix  notation, 

<R.R|>  =  4  R2(s..-^)«RR>=1  R2a-SS)] 


Thus  we  obtain 


/d2f  \(R)  -  IM  -  d£ 

W  /secular"  R  “  dt  ' 

Now,  if  the  satellite  equation  ignoring  relativity  was 


df£ 

dt2 


=  F 


(5-14) 

(5-15) 

(5-16) 

(5-17) 

(5-18) 


then  GRT  changes  it  to 


—  zn  x  =  f 
dt2  dt 


def  3M  - 
=  2R  w 


(5-19) 

(5-20) 


If  we  now  transform  Eq.  (5-19)  to  a  system  of  axes  centered  in  the  earth  and  rotating  with 
—  2 

angular  velocity  ft,  then,  neglecting  terms  proportional  to  ft  ,  one  obtains,  in  terms  of  the 
moving  axes,  just  Eq.(5-18).  So  the  De  Sitter-Schouten  Effect  is  that  of  a  secular  rotation  of 
the  nonrelativistic  orbit  of  a  satellite  by  an  amount  which  depends  only  on  the  distance  of  the 
primary  to  the  sun.  For  the  case  of  the  moon,  this  effect  amounts  to  about  1.9  seconds  of  arc 
per  century.  It  therefore  may  very  well  be  observable. 
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VI.  PROPAGATION  OF  LIGHT  RAYS  IN  THE  SOLAR  SYSTEM 


We  will  make  the  usual  assumption  that  the  trajectory  of  a  light-ray  is  the  same  as  that  of  a 

massless  particle  (photon)  so  that  the  equations  of  motion  are  the  equations  of  a  geodesic  corre- 

4 

sponding  to  a  massless  particle  (null  geodesic): 


,2  i  .  ,  k  ,  ft 

d  x  ,  r  1  dx  dx  _ 

2  ki  da  da 

da 


,  i  ,  k 
dx  dx 

^ik  da  da 


[Eq.  (2-17)] 


These  equations  will  determine  both  the  path  and  the  parameter  a  which  is  called  a  "distin¬ 
guished  parameter  along  the  geodesic."  Since  the  concept  of  a  distinguished  parameter  is  impor¬ 
tant  in  our  discussion  of  the  Doppler  Effect  (see  Sec.  VII),  we  will  expand  somewhat  on  the  origin 
of  the  last  equations  and  on  geodesics  in  general. 

A  non-null  geodesic  can  be  defined  as  a  curve  along  which 

6  CB  ds  =  0  [ Eq.  (2-11)] 

JA 


where 


(ds)  =  — g.  .dxldx^ 
If  we  parameterize  the  curve  by 
x*  =  x*(a) 

then  we  want  to  make 


[Eq.  (2-12)] 


r 


y  dx  dxJ 
*ij  da  da 


da  =  0 


(6-1) 


(6-2) 


The  Euler- Lagrange  equations  for  this  problem  lead  to 

k  ...l 
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dZxL  Fi  dx^  dx*  _  _d_/in  dx* 

^2  jk  da  da  daV  da  da 


assuming  that 


is  ^  ° 

da 


(non-null  condition) 


(6-3) 


(6-4) 


Since  the  parameter  a  is  arbitrary,  we  can  choose  it  in  such  a  way  as  to  simplify  the  equa¬ 
tions.  We  choose  it  to  be  proportional  to  the  arc  length  along  the  non-null  geodesic  by  setting 


dxL  dxJ  _  ± a2 
gij  daD  d<rD 

where  a  is  an  arbitrary  positive  number,  in  which  case  Eqs.  (6-3)  reduce  to 

d2x*  „  i  dx*  dx-*  _  q 
dff2  jk  daD  daD  ‘ 


(6-5) 


(6-6) 


Any  cD  which  leads  to  Eq.  (6-6)  is  called  a  distinguished  or  special  parameter  along  the  geodesic. 
We  see  immediately  that  if  a^  is  distinguished,  so  is 
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°D2  =  AcrDl  f  B 


(6-7) 


where  A  and  B  are  arbitrary  constants.  In  particular,  we  can  always  choose  a^  so  that 


dx*  dx^  _  ±1 
gij  daD  daD  " 


(6-8) 


Except  for  the  sign,  this  looks  just  like 

ds2  =  —  g.jdx^dx^  .  [Eq.(2-9)] 

The  two-fold  sign  appears  in  Eq.(6-8)  but  not  in  Eq.(2-9)  because  Eq.(2-9),  though  traditional, 

is  ambiguous.  The  right-hand  side  of  Eq.(2-9)  can  be  either  positive  or  negative  whereas  the 

arc-length  is  conventionally  a  real  number.  For  a  curve  which  represents  the  path  of  a  nonzero 

mass  particle,  the  right-hand  side  of  Eq.  (2-9)  is  positive  so  there  is  no  problem  there,  and  it  is 
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tacitly  understood  that  in  other  cases  an  extra  minus  sign  will  be  inserted  when  necessary  for 
consistency. 

The  above  analysis  breaks  down  for  null  curves,  i.e.,  curves  along  which 

ds2  =  — g^dx^dx^  =  0  (6-9) 


since  the  right-hand  side  of  Eq.  (6-3)  is  then  meaningless.  In  order  to  modify  the  formalism  to 
handle  this  case,  we  replace  the  variational  principle  of  Eq.  (6-2)  by 


*  C  dx  dx'1  ,  n 

6  \  g.  -t - da  =  0 

J  6i]  da  da 


(6-10) 


where  a  is  no  longer  completely  arbitrary  but  is  restricted  by  the  subsidiary  condition 


dx*  dx^ 
^ij  do-  da 


(6-11) 


This  variational  principle  leads  to  Eq.  (6-6)  and,  because  of  Eq.  (6-11),  o  can  be  identified  with 
the  arc-length  when  the  geodesics  are  non-null  curves.  These  curves  which  satisfy  Eq.  (6-6) 
and  which  are  null  curves  are  called  null  geodesics.  One  can  show  that  Eq.  (6-6)  implies  that 


A.L  ,<i>L^h  =  o  (6-12) 

da \6ij  da  da  / 

so  that  the  expression  in  Eq.  (6-11)  has  a  constant  value  along  a  given  geodesic.  As  long  as  this 
number  is  nonzero  we  can  normalize  its  absolute  value  to  unity  and  then  a  becomes  the  arc- 
length  along  the  curve.  It  also  follows  from  Eq.  (6-12)  that  a  geodesic  which  is  null  anywhere 
is  null  everywhere,  a  can  no  longer  be  identified  with  the  arc-length  along  a  null  geodesic.  We 
have  thus  derived  Eq.  (2-17).  If  we  wish  to  use  a  nondistinguished  parameter  along  a  null  geo¬ 
desic  then  we  replace  Eq.(2-17)  by  the  more  general  form 


d2x*  pi  dx‘i  dx^  _  d2\  /  cUv2  dx* 
dx2  +rjk  dX  dX  '  d(JV  d"  ^ 


(6-13) 


g 


dx^ 
ik  dX 


(6-14) 
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(6-15) 


In  particular,  if  we  choose  the  time  as  the  parameter,  Eqs.(6-13)  and  (6-14)  become 

(or  =  1,  2,  3) 
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d  x 
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dx 
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dx 
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dt 

dt 
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a  1 
q 

ro| 

dx* 

dx^ 

gik 

dt 

dt 

=  0 

G  dtv2  dx^ 
da  da 


(6-16) 


(6-17) 


It  can  be  shown  that  Eq.  (6-16)  is  redundant,  i.e.,  if  Eq.  (6-15)  and  Eq.  (6-17)  are  satisfied  then 
Eq.  (6-16)  holds  automatically.  Thus  we  need  only  consider  Eqs.(6-15)  and  (6-17). 

The  geodesic  Eqs.(6-15)  and  (6-17)  involve  T  and  g^.  The  g.  .,  up  to  the  post-Newtonian 
terms,  are  given  in  Eqs.  (2-68)  to  (2-70).  It  will  turn  out  that  we  can  ignore  the  post-Newtonian 
terms  as  their  effect  is  so  small  as  to  be  undetectable  at  present  and  in  the  near  future.  Thus, 
we  can  write 


gcv4  ~  ° 


(6-18) 

(6-19) 


g44  =  -  1  -  ^ 


(6-20) 


where  we  have  used  units  such  that  c  =  G  =  1  and  $  is  twice  the  Newtonian  potential 


ip(x) 


def  ^ 

=  -2  E 


|X-X.| 


(6-21) 


and  where  we  have  taken  the  origin  of  our  coordinate  system  in  the  sun  since  to  this  order  the 
velocity  of  the  sun  can  be  ignored. 

The  lowest  approximation  is  to  ignore  the  gravitational  potential  altogether,  i.e., 


>>  =  +6 

ga/3  a/3 


*4?  =_1 


other  g. .  =  0 


This  causes  all  the  t\  to  vanish  and  leads  to 
JK 


(6-22) 


d2X°' 

dt2 


d2t  /  _ctt\2  dx_ 

,  ZKd(j)  dt 
da 


a 


(6-23) 


dx**  dx**  _  q  (6-24) 

dt  dt 

Differentiating  Eq.  (6-24)  and  combining  with  Eq.  (6-23),  we  find 

d2t  .  0  (6-25) 

da 

and 

d2x**  _  n  (6-26) 

2 

dt 
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so  in  this  approximation  the  distinguished  parameters  and  geodesics  are 
cr  =  At  +  B 

and 


(6-27) 


x*  =  Dat  +  Ea 


(6-28) 


where  A,  B,  and  Ea  are  constants.  Light  travels  in  a  straight  line  and  since  in  our  units 
the  speed  of  light  far  away  from  matter  is  unity,  we  write 


xa(t)  =  x"(t  )  +  (t  —  tQ)  n01 
o 


(6-29) 


where  na  are  the  components  of  a  unit  constant  vector 


a  a  a  /\  . 

n  n  =  n  •  n  =  1 


(6-30) 


as  the  first  approximation  to  the  path  of  a  light  ray  which  passes  through  the  point  with  coordi¬ 
nates  xa(t  )  at  time  tQ  and  is  traveling  in  the  direction  of  the  unit  vector  n.  We  can  ignore  the 
distinction  between  covariant  and  contravariant  indices  on  space  components  of  vectors  since 
the  metric  is  Minkowskian  to  this  order,  i.e.,  is  given  by  Eq.  (6-22). 

We  now  find  the  corrections  to  this  trajectory  in  the  next  approximation 


xa  (t)  =  x^(t)  +  xa  (t) 
1  o 


(6-31) 


The  Christoffel  symbols,  T  ?  ,  are  to  first  order  in  ip, 

JK 
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1  _  i_  _d£_ 

Ox3 


ri4  =  °' 


etc. 


Using  Eqs.  (6-32)  and  (6-29)  for  x(t),  we  can  convert  Eqs.  (6-15)  and  (6-17)  into 

o 


(6-32) 


-^5  xa  =  2  Jt  A®  -  ii- 

dt2  i  7? 


dx 


(6-33) 


and 


a  d  a 
n  dt  *  " 


(6-34) 


In  vector  notation,  these  become 

,2 

-S  X  =  2n(n  •  V)#-V4> 
dt  1 


(6-35) 


and 


n  •  dt  * =  * 


(6-36) 


34 


Fig.  1.  Special  coordinate  system. 


These  equations  can  be  integrated  directly.  However,  to  make  things  more  transparent,  let 
us  look  at  a  ray  in  the  x-y  plane  which  in  lowest  order  moves  in  the  positive  x-direction.  This 
obviously  involves  no  loss  of  generality  and,  as  we  will  see,  our  results  can  be  transcribed 
easily  into  coordinate-free  vectorial  notation.  Referring  to  Fig.  1,  let  the  ray  start  out  at 
(x(t),  y ( r ) )  at  time  r  with  n  =  (1,  0).  Then 


x(t)  =  x( r )  +  (t  -  tq) 
o  o 


x(t)  =  1 
o 


y(t)  =  y(r) 
o  o 

y(t)  =  0  . 

o 


The  equations  of  motion  for  the  next  approximation  become 

= s 

1 

y(t)  dip 

1  ay 

x(t)  =  ip 
1 

Integrating  the  equations 


<,T  . 


y(t)  .-f* 

1  “T 

y(t)  =  -  C  dT  \  1  dT  +  y(r)  (t  -  r)  +  y(r) 

1  ^  1  1 

=  -  j  (t  -  T)  dT  +  y(r)  (t  -  t)  +  y^r) 

?(t)  =  f  4>( T)  dT  +  x(t) 

JT  i 

where  the  integrals  are  evaluated  along  the  zeroth-order  trajectory,  i.e.,  if 
F  =  F(x,  y,  z,t) 

\  F(T)  dT  d=f  (  F  I  x(T),  y(T),  z(T),  T  dT  . 

d  d  l()  O  O  > 


(6-37) 

(6-38) 

(6-39) 

(6-40) 

(6-41) 

(6-42) 


(6-43) 


(6-44) 


(6-45) 
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Adding  on  the  zeroth-order  terms,  we  obtain  the  basic  formulas 


x(t)  =  x(t)  +  x(t)  +  (t  —  t)  +  \  4i T)  dT 
o  1 

y(t)  =  y(r)  +  y(r)  +  y(T)  (t  t)  1  (t  -  T)  (T)  dT 
oil  JT  9y 

x(t)  =  1  +  t) 

y(t)  =  y (t)  (T)  dT 


(6-46) 

(6-47) 

(6-48) 

(6-49) 


We  now  evaluate  the  integrals  in  the  above  equations  for  the  part  of  4>  which  is  due  to  the  sun 
since  we  will  show  later  that  this  is  all  we  need.  There  is  no  loss  in  generality  in  placing  the 
sun  at  the  origin  of  our  coordinate  system  so  that 

4(x,y,z)  =  —  -  ZM  (6-50) 

/  2  ^  2  2 

n lx  +  y  +  z 

and 


We  find 


4>(  T) 


-2M 


J 


(x(r)  +  (T  -  r)]2 
o 


+  [y(r)l2 
o 


(6-51) 


(H T)  dT  =  — 2M  \ 


dT 


r  { ( x( t )  +  (T  -r))2  +  y2(r)}1/2 


—  2M 


x(t ) 

r°  _ do _ 

^x(t)  [O'  2  +  y2(T)]1//2 


=  -2M  In 


(X2(t)  +  y2(r)]l/2  +  x(t) 


(x2(r)  +  y2(r)l1/2  +  x(t) 
o  o 


(6-52) 


r 


|y  dT  =  2M  y(r)  \ 


x(t) 

nO 


dev 


2M 


\  ~t  2  7  2,  m3/2  y(T)  r  277%  2,  ..1/2 

O  jx(t)  [o  +  y  ( T )]  ^  [x  (t)  +  y  (r)] 


x(t) 

o 


x(t) 

o 


.  2.  .  2.  ..1/2 

[x  (r)  +  y  (r)] 
o  o 


(6-53) 


T 

£ dT.  iT  1 m  dT  = 


2M 
y  <t) 
o 


[X2(t)  +  y2(r)]l/2 


(x2(r)  +  y2(r)]l/2  + 


-(t-r)x(r)  1 

o 

TJl 


[x2(t)  +  y2(T>] 
o  o 


(6-54) 
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Therefore  we  finally  obtain 


and 


x(t)  =  x(t)  +  (t  —  t)  —  2M  In 
o 


[X2(t)  +  y2(r)ll/2  +  x(t) 


[x2(r)  +  y2(r)]l/2  +  x(r) 
On  o 


+  x(t) 
1 


y(t)  =  y(T)  +  (t  -  t) 


y(r)  +  2.  ,  .  2.  ~T7 2 

1  J  [X  (t)  +  y  (t)] 


2M 


x(t) 
o 


2M  /r  +  \  x  2,  m1/2 

“  TTTj  ([x  (t)  +  y  (t)] 

7  o  o 

o 


x(t )  =  1  — 


r  2.  »  2.  v,l/2, 

-  [x  (t)  +  y  (t)]  } 

o  o 

_ 2M 

[X2(t)  +  y2(x)] 1/2 


y(t)  + 


x(t) 

O 


x(  T  ) 

O 


2M  o  -  \  x  2M  o 

y(r)  .  2,t.  ,  2,  ..l/2  "  y'T'  +  y(r)  .  2,  .  .  2,  ..1/2 

£  [x  (t)  +  y  (t)]  o  fX  T  +  y  T" 


(6-55) 


(6-56) 

(6-57) 


(6-58) 


We  can  make  the  content  of  these  formulas  clearer  if  we  observe  from  Fig.  1  that  y<T>  is  just  "d," 

the  distance  of  the  closest  approach  of  the  undeflected  ray  to  the  sun,  and  that[xQ(t)  +  y  ( r )]  u 
is  just  the  distance  of  the  point  on  the  undeflected  path  from  the  sun: 


r(t)  +  x(t) 

x(t)  =  x( r )  +  (t  -  t)  -  2M  In  r(T)  +  x(~)  + 


(6-59) 


y(t)  =  d  +  (t  —  t) 


x(t) 

•  /  v  ,  2M  o _ 

y(T)  +  —  7(7) 


—  [r(t)  -  r(r)] 


x(t)  =  1 


2M 

r(t) 


x(t)  x(t) 

+  \  X  2M  °  •/  N  ,  2M  o 

y(t)  +  ~  7(tj  =  y(T)  +  “d"  7(7)  • 

These  formulas  will  now  be  applied  to  obtain  the  results  of  interest  to  us: 


(6-60) 

(6-61) 


(6-62) 


A.  APPARENT  SPEED  AND  ABERRATION  OF  LIGHT 

For  an  observer  at  rest  relative  to  the  sun,  the  apparent  speed  of  light  c  and  the  angle  a 
which  the  ray  will  make  with  the  x-axis  are  obviously  determined  by  the  equations 


c  cos  a  =  x  =  1  +  x 
1 


c  sin  a?  =  y  =  y 

1 


=  /( 1  +  x)2  +  y 

\ I  1  1 


•2  4  x  * 

y  »  1  +  x 

1 


(6-63) 

(6-64) 

(6-65) 
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(6-66) 


tana  =  y/(i  +  x)  ~  y  ~  cv 
1  1  1 


B.  DEFLECTION  OF  STARLIGHT  SEEN  FROM  THE  EARTH 

Suppose  the  light  ray  arrives  at  the  earth  at  the  time  t: 


X<T)  =  xe(t) 

o 


x(t)  =  0 
1 


y  (t)  =  d 

o 


y(T)  =  0 
1 


(6-67) 


Furthermore,  let  the  light  appear  to  come  along  the  x-axis,  i.e., 

y(r)  =  0 
1 


(6-68) 


Then  the  actual  velocity  of  the  light  at  any  time  is,  from  Eqs.  (6-61)  and  (6-62), 

<  2M 

x(t)  =  1  ~  rftj 

,  2M  o_]_  _  2M  xE(t) 
d  r(t)  '  d  rE(r) 


(6-69) 


(6-70) 


and,  by  Eq.  (6-66), 


o(t)  = 


x  (r)  x(t) 

2M  XE m  _  2M  o _ 

d  r  jp(  t  )  d  r(t) 


As  t  -►  —  x(t)  —  oo  and  r(t)  -►  |x(t)|  so 
o  o 


(6-71) 


2M  .  2M  xE(t) 

a(-oo)  =  lim  »(t)  =  -7-  +  -T-  — j— v 

t— -  d  d  rE  T 


(6-72) 


|  3  63-5619  I 


APPARENT  PATH 


Fig. 2.  Path  of  light  ray. 


We  illustrate  this  result  pictorially  in  Fig.  2  for  the  case  where  the  earth  and  the  star  are  on 
opposite  sides  of  the  sun.  Note  that  the  deflection  increases  monotonically  along  the  path  and 
that  the  deflection  between  t  =  —  °°  and  the  closest  approach  to  the  sun  is  2M/d.  The  ray  must 
approach  the  sun  closely  for  the  bending  to  be  significant  since 


M 

R 


0  V  45 


(6-73) 


So  if  we  assume  x^  »  d, 


I  \  4M  lR  @\ 

„(-«)  „  —  =  1 V 7 5  {^) 


(6-74) 


which  is  a  classic  result. 
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c.  ABERRATION  OF  STARLIGHT  RELATIVE  TO  ABERRATION  OF  A  PLANET 
SEEN  FROM  THE  EARTH 


Suppose  a  planet  and  a  star  appear  to  be  along  a  straight  line,  i.e.,  their  images  coincide  as 
seen  from  the  earth.  If  the  light  rays  have  passed  near  the  sun  on  their  way  to  the  earth  then 
their  directions  have  changed  relative  to  what  they  would  have  been  had  they  traveled  in  straight 
lines.  We  now  calculate  this  shift. 

Let  us  first  assume  that  the  earth  is  at  rest  relative  to  the  sun  and  again  let  the  light  appear 
to  come  along  the  x-axis.  Then  we  can  use  Eq.  (6-72)  for  the  starlight: 


2M  2M  xE(t) 

01  st  d  d  rE(r) 


[Eq.  (6-72)1 


Let  the  light  leave  the  planet  at  time  r  —  t.  .  If  it  appears  to  come  along  the  x-axis,  then  by 
Eqs.(6-61)  and  (6-62) 


.  2M  xp(T  Tl> 

y(r-Ti)  +  —  ^ 


3(r  -Tl> 


2M  xE(t) 
d  rE(r) 


(6-75) 


,  .  2M 

“p*7"7!*  =  —  r 


xE(t)  _  2M  Xp(T 


Jr) 


71)  +  2M 
rp(T"Tl)  rp(T_Tl) 


(6-76) 


The  last  term  can  obviously  be  dropped  since  it  is  ~d/r  as  big  as  the  other  two  and  so  com¬ 
bining  Eq.  (6-76)  and  Eq.  (6-72)  we  obtain 


2M 

ast  %  "  d 


xn(r  “Tl) 

1  +  — k-? - ^ 

rD(T  "Tl) 


(6-77) 


This  is  not  a  very  practical  form  for  the  result,  since  we  would  like  to  express  it  directly  in 
terms  of  the  relative  coordinates  of  the  planet  and  of  the  earth  rather  than  in  terms  of  "x  "  and 
"d."  Consider  the  straight  line  passing  through  the  points  ?1  and 


def 


'(\)  =  +  X(?2  -  ?1)  =  +  XR 


-00  <  X  <  °o 


?(0)  =  F 


1 


F(D  =  ?z 


The  distance  A(X)  between  the  point  F(X)  on  the  line  and  the  sun  (at  the  origin)  is  given  by 

A2(X)  =  |  rr1  +  XR  |2  =  r^  +  2xR  •  +  X^R2 

The  distance  of  closest  approach  is  obtained  by  differentiating  Eq.  (6-80): 
k2  _ 


0  = 


dA 

dX 


2R  •  +  2XR 


at  minimum 


R  •  r . 


R 


Now 


2  2 


d2  =  A2  .  |(r.  X  r ?)\z  = 

mm  R2  1  2  R 


rl  r2  2 

- —  sin  4  (r1,  r^) 


(6-78) 

(6-79) 

(6-80) 


(6-81) 


(6-82) 
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and 


where 


d  =  |  sin  4  (r^,  r2)  | 


R  =  |r1  -  r2| 


(6-83) 


We  now  can  determine  the  quantities  "x"  and  Md”  for  a  straight  line  passing  through  the  positions 
of  the  earth  at  time  r  and  of  the  planet  at  time  t  —  t ^  in  terms  of  the  coordinates  of  the  earth  and 
of  the  planet  in  the  earth-sun-planet  plane: 


ri?  ~  cosG.,  r^  sin  9^)  ;  r  =  (r  cos  0  ,r  sin0  )  (6-84) 

xl.  IL,  HilL  £L  pp  p  p  p 

rZ  =  rE  +  rp  “  2rErp  cos  (0E  ~  6p)  (6‘85) 


d 


E 

R 


|  sin(0E  -  0p)  | 


x 


P 


r 

P 


=  ^  [rE  cos(eE-ep)-rp] 


(6-86) 

(6-87) 


and  the  angle  a  is  measured  counterclockwise  from  R. 

In  the  above  analysis  of  the  deflection  of  light,  we  have  assumed  the  earth  to  be  at  rest 
relative  to  the  sun.  When  we  take  the  earth's  motion  into  account,  we  find  that  from  the  view¬ 
point  of  an  observer  on  the  moving  earth  the  direction  a  ray  of  light  appears  to  be  traveling 
differs  from  the  direction  that  would  be  seen  by  a  stationary  observer  at  the  same  point  in  space. 

This  effect  is  called  aberration  and  is  calculated  by  making  a  Lorentz  transformation  from  the 
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sun's  rest  frame  to  the  instantaneous  rest  frame  of  the  observation  point  on  the  earth.  Since 

-4 

the  relative  earth-sun  speed  is  about  10  c,  and  angle  of  arrival  measurements  are  limited  in 

precision,  it  is  sufficiently  accurate  to  use  the  lowest  order  aberration  correction,  which  may 
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be  expressed  as  follow's: 


\  -  Ae  +  (Ae  •  VE)  VE  dtf  Aq  (6-87a) 

w'here  A  and  A  are  vectors  in  the  direction  of  the  ray  in  the  rest  frames  of  the  earth  and  sun 
e  °  -► 

respectively,  and  is  the  velocity  of  the  earth  relative  to  the  sun  (the  fact  that  Ap  and  Aq 
have  different  magnitudes  is  irrelevant  since  we  are  only  interested  in  the  directions  along 
which  they  point).  Wrhen  Eq.  (6-87a)  is  combined  with  the  gravitational  deflection  wre  obtain 


A  /A  A\  A 

-  /\  2M  rE  (rE  '  e)  6  .  ,/N  a  ,  r, 

e  =  e  +  -  -  +  (e  •  V„)  Vt 

o  r^  .a  a  E  l 

E  1  —  r^  .  e 


( 6-87b) 


e  unit  vector  along  direction  of  ray  in  earth's  rest  system 


-►  def 

eQ  =  vector  along  direction  of  ray  in  sun's  rest  system 
d  cf 

?E  =  vector  from  sun  to  earth 
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We  can  rewrite  Eq.(6-87b)  for  a  geocentric  coordinate  system.  Let 

n  d=f  unit  vector  from  earth  to  point  in  sky  from  which  light 
appears  to  be  coming 


def 

r  =  vector  from  earth  to  sun 
s  -  — 


def 


velocity  of  sun  relative  to  earth 


( 6-87c) 


Then  we  make  the  following  substitutions  in  Eq.  (6-87b) 


A  A 

e  —  n 


r„  -  r 
E  s 


( 6-8  7d^ 


V  —  —  V 

V  E  Vs 


and  arrive  at  the  final  formula 


A  .A  A >  A 

r  —  (r  •  n)  n  ~ 

— ►  A  .  2M  S  S  /A  -rr  V  Tr 

n  =  n  + - t  (n  •  V  )  V 

o  r  a  a  s  s 

s  1  —  r  •  n 
s 


( 6-87e) 


D.  TIME  OF  FLIGHT  OF  LIGHT 

Suppose  we  have  two  bodies,  say  the  earth  and  a  planet,  whose  trajectories  FE<t)  and  Fp(t) 
are  assumed  known.  We  send  light  (or  a  radar  beam)  from  F£  at  a  time  t1  and  want  to  know  at 
what  time  t1  +  T  it  will  reach  the  planet.  If  GRT  could  be  ignored,  we  would  aim  the  beam  along 
the  straight  line  from  F^t^  to  F  (tj  +  T)  and  the  light  would  travel  along  this  line  with  a  speed 
of  unity  (in  our  units).  Thus  the  trajectory  of  the  beam  in  flat  space  would  be 

F  (t,  +  T)  -  rF(t,) 

F(t)  =  ?rp[tA)  +  — -  (t  -t.)  (flat  space)  (6-88) 

Irp(t1  +  T)-?E<t)| 

and  the  time  of  flight,  T,  must  be  obtained  from  the  implicit  equation 


rp^l  +  "  rE^tl ^  + 


rytj+T)-?^) 

I  rp^l  +  -  rE*tl^ 


T 


(6-89) 


or 


T  =  l?p(t1  +  T)  -  r’E(t1)|  (flat  space) 


(6-90) 


GRT  introduces  two  corrections:  The  actual  path  is  not  a  straight  line  and  the  local  speed 
of  light  is  not  unity.  This  can  be  seen  from  Eqs.(6-59)  to  (6-62)  which  automatically  take  these 
effects  into  account.  Choosing  the  x-axis  to  connect  Fgft^)  to  F  (t^  +  T), 


x(t)  =  Xgft^  +  (t  —  t^)  — 2M  In 


r(t)  +  x(t)  1 

rE(t1)  +  xE(t1)J 


(6-91) 


and  the  time  of  flight,  T,  is  to  be  determined  from 


x  (t,  +  T)  =  Xp( t , )  +  T  — 2M  In 

pi  XL  1 


rp(ti  +  T>  +  xP(ti  +  T) 

rTr,(t . )  +  x^(t . ) 


lEur 


E  1 


(6-92) 


41 


The  last  equation  can  be  rewritten  in  terms  of 


— ►  def  — 

^a'V  =  VV  "  rE(tb} 

n(ta,  tfa)  d=f  R/|R| 


(6-93) 

(6-94) 


xp(tl  +  T)  = 


R(tj  +  T,  td) 


R 


V*!  +  T) 


(6-95) 


xE<tl )  ~ 


R(tt  +  T,  tt) 
R 


rE(tl) 


(6-96) 


as  (Ref.  65) 


T  =  R(t,  +  T,  t,)  +  2M  In 
1  1 


yti +  t>  + « •  y*i +  t> 


rE(V  +  n  ’  ?E(tl} 


(6-97) 


which  is  to  be  contrasted  with  the  flat-space  equation  (6-90). 

If  the  signal  is  reflected  back  from  the  planet  so  that  it  reaches  the  earth  at  time  t^  +  T  +  T^, 
then  Eq.  (6-59)  gives  for  the  return  trip 


x(t)  =  x  (t,  +  T)  +  (t- T  -  tJ  -2M  In 
pi  l 


r(t)  +  x(t) 


rp<t1  +  T)  +  xp(tl  +  T) 


(6-98) 


where  the  positive  x-axis  is  now  in  the  direction  of  r^ft^  +  T  +  T^)  —  r  (t^  +  T),  i.e.,  essentially 


reversed: 


+  T  +  T . )  —  r  (t.  +  T) 

x(t)  =  - - - - -  •  ?U)  =  — n(t .  +  T,  t.  +  T  +  T.)  •  r(t) 

|rE(t1  +  T  +  T4) -?p(t1  +  T)| 


(6-99) 


so 


Tj  =  R(t1  +  T,  +  T  +  Tt)  +  2M  In 


r£(t1  +  T  +  Tt)  -n*  •  +  T  +  Tj) 


rp(tl  +  T)  — n1  ■  ?p(tt  +  T) 


(6-100) 


def 


n'  =  n(t .  +  T,  t,  +  T  +  T.) 

11  1 


(6-101) 


Note  that  the  arguments  of  the  logarithms  in  Eqs.  (6-97)  and  (6-100)  are  essentially  equal  as  a 
consequence  of  the  identities 

(r  +  n  •  ?  )  (r  —  n  •  7  )  =  r  2  —  (n  •  7  )2  =  (7  x  n)  •  (r  x  n) 

P  PP  PP  PP  P 

(7P  X  7E)Z  j2 
=  -  =  d 


and 


(rE  +  n  .  rrE)  (r£  -  n  •  r^,)  =  rE  -  (n  •  r£)2  =  (?E  x  n)  •  (r^,  X  n) 


(fp  *  ?E)2  .2 


=  d 


R 


(6-102) 
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The  GRT  one-way  time  delay  is  thus  approximately  doubled  for  a  round  trip.  In  terms  of  the 
distance  of  closest  approach  to  the  sun,  the  outbound  GRT  delay  is 


(fit)  =  2M  In 
out 


r  ± 
P 


/ 7*-?] 

v  p 


-J*E~ 


(6-103) 


Use  of  plus  or  minus  is  determined  by  whether  the  earth  and  planet  are  on  the  opposite  or  same 
side  of  the  sun,  respectively.  This  is  obviously  greatest  near  superior  conjunction  and  then  we 
can  write 


(^)  • 


(6t)  .  «  2M  In 

out 


(6-104) 
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For  d  =  R  ,  this  term  is  about  10  sec  for  Venus  and  Mercury.  We  can  also  rewrite  (fit)  ,  in 
o'  J  out 

terms  of  the  earth-planet  coordinates  in  the  earth-planet-sun  plane  as  given  in  Eqs.  (6-84)  to 
(6-86): 


(fit)  ,  =  2M  In 
out 


r^  +  r  +  R 
E  p 


r„  +  r  -  ft 

E  p 


R  =  (rE  +  rp2-2rErP  cos(eE-ep)]l/2 


(6-105) 


The  expression  for  the  return  trip  has  exactly  the  same  form  but  the  various  quantities  are 
evaluated  at  different  times.  For  later  use,  we  also  calculate: 


-£-(fit)  .  =  2M 

dt  out 


r_  +  r  +  R  r^  +  r-R 
E  p  _  E  p 

r-,-.  +  r  +  R  r  ^  +  r  —  ft 


E  '  Ap  ’  iV  E  "p 
which  for  the  case  of  coplanar  orbits  becomes 

2M 


( 6- 1 06a) 


37  ( fit)  ,  =  33-  ( r'  r  —  r  r^) 
dt  out  R  E  p  p  E 


(^lp) 

\  rprE  / 


+  -3-(rE  +  rp)  (®E  _  ®p)  tantr<eE  “  ep)1  •  (6-106b) 


<  - 1 2 

The  first  term  is  always  small  ( — 10  )  but  the  second  term  becomes  much  larger  near  superior 

conjunction.  We  then  have 


.“'iIOE  —  9p|l“t  (r^Y*  ■  ‘‘'“e-9, 

P 


P  \l  £ 


(>)* 


(6-107) 


and 


(6t)  .  «  ±4M  - 

dt  out  rn  +  r 


r^r  0—0 
E  p  E  p 


(6-108) 


which  has  the  values,  for  coplanar  orbits 


4(6t)  .  »  T2.1  X  1 

dt  out 


=f6.9  x  1 


•-(?) 


Venus 


Mercury 


(6-109) 
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The  basic  equations  (6-97)  and  (6-98)  are  solved  by  an  iteration  scheme  using  the  best  known 
values  of  r^ft)  and  r^t).  Setting  t^  =  0,  which  entails  no  loss  of  generality,  the  equations  to 
solve  are 


T  =  F(T,  0) 
T'  =  F(T,  T') 


F(tA-  tB)  "  R(tA’  tB)  +  A(tA*  tB) 

A’ 


R(tA’  tB>  lrp^A^  rE(tB) 


A(tA,  tB)  =  2M  In 


rE(tB)  +  rp(tA)  +  R(tA’  tB) 


R(t a.  tn) 


P'  A 


A’  B 


We  then  iterate 


(6-110) 


(6-111) 


Tt  =  F(0,  0)  • 

t2  =  F(T4,  0)  , 
t3  =  F(T2,  0) 

=  f(t3,  t3) 

t2  =  f(t3,  t3  +  T\)  ,  . 

t3  =  f(t3,  t3  +  t2) 

To  show  that  this  is  sufficiently  accurate,  let  us  look  at  Eq.  (6-112)  and  Eq.  (6-110), 
T  -  tn+i  =  F(T*  0)  “  F(TN'  0) 


(6-112) 


(6-113) 


Then 


OF 

at  A 


N  '  B" 


def 

-m  (T-Tn)  =  (T-T„>  a 


N' 


(6-114) 


where  0  ^  a  <  1  and  where  we  have  used  the  mean  value  theorem  of  elementary  calculus.  This 
can  be  rewritten  as 


T  -  T 


N+l 


°<tn-h  tn) 

1  -  a 


(6-115) 


By  repeatedly  applying  the  mean  value  theorem  and  using  the  fact  that  the  time  derivative  has 
the  same  order  of  magnitude  throughout,  we  conclude  that 


T  - 


T 

1  N+l 


N+l 
a _ 

1  —  a 


(6-116) 


Now,  T 


1 


1 000  sec  and 
dR 


R 


dt 


=  (V  -V  )  . 

K  p  V  E  R 


(6-117) 


so  if  we  combine  Eq.  (6-117)  and  Eq.  (6-106)  we  find  2  X  10~4  and  |  T  —  T.J  <  4  X  10~9  sec. 
The  time  of  flight  is  thus  given  by  T^  +  T^  to  within  10~  sec.  In  view  of  Eq.  (6-106)  it  is  not 
really  necessary  to  iterate  A  more  than  once;  ^(T^  0)  can  be  used  throughout. 
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We  have  ignored  the  effects  of  the  gravitational  fields  of  the  planets  on  the  time-delay.  Let 
us  show  that  these  terms  are  in  fact  negligible.  Going  back  to  Eq.  (6-46)  and  noting  that  the 
closest  we  come  to  the  center  of  a  planet  is  one  planetary  radius  away,  we  see  that 


(<5t)  -  .  ,  «4M  In 

planet-earth  p 

round  trip 


+  4Me  In 


(6-118) 


-9 

which  is  a  time  delay  of  less  than  10  seconds  for  Venus  and  Mercury. 

Finally,  we  wish  to  relate  the  value  of  the  time-of-flight  we  have  obtained  to  that  which 
would  be  measured  by  an  observer  on  the  earth  in  terms  of  an  atomic  clock  at  the  radar  site. 
The  invariant  interval  between  two  events  is  given  by  Eq.  (2-9): 


ds2  =  — g^dx^dx^  [Eq  (2-9)] 

which,  in  our  approximation,  becomes  [cf.  Eqs.(6-18)  to  (6-21)]: 

ds2  =  (1  +  ip)  dt2  —  (1  —  i/>)  (dx2  +  dy"-  +  dz2) 


=  dt2  (i  r<^)2  +  <^)2  +  <^)2i) 

7  o  def  ? 

~  dt  {l  +  ip  —  v  }  =  dT  .  (6-119) 

A  fundamental  postulate  of  GRT  identifies  this  quantity  with  the  square  of  the  proper-time  interval, 
i.e.,  the  time  measured  by  a  clock  at  rest  relative  to  the  moving  body.  Put  another  way,  if  a 
clock  at  rest  and  very  far  away  from  all  other  masses  would  read  "dt"  then  an  identical  clock  on 
the  earth  would  read  "dT."^  The  atomic  clock  on  the  earth  then  would  measure  a  ”proper-time- 
of-f light”  of 

1 1  +Tf 

r  =  f  [1  +  4>  -  v2]l/2  dt  (6-120) 


if  tj  is  the  coordinate  time  at  which  the  signal  leaves  the  earth. 

The  integrand  in  Eq.  (6-120)  is  not  constant  due  to  the  motion  of  the  various  bodies  in  the 
solar  system  and  to  the  rotation  of  the  radar  site  with  the  earth,  so  the  relations  between  the 
coordinate-time  and  the  local  atomic-time  epochs  are  quite  complicated.  Since  M.  Ash  and 
I.  Shapiro  have  studied  this  point  extensively,  we  will  not  discuss  it  here  but  will  merely  show 
that  the  variation  of  the  integrand  during  the  time-of-flight  can  be  ignored,  i.e.,  we  can  use 

T  =  (l  +  0(t1)  —  [v(t1)]2}  Tf  .  (6-121) 

The  next  term  in  a  Taylor  Series  expansion  of  t  would  be 


dij) 

dt 


Tr  -  v-^f 


[1  + 


dv 

dt 

J7z 


Now, 
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(6-122) 


r  =  r  +  p 
c  r 


v  =  \rc  +  c j  X  p 


a-a^+ajX(5xp) 


where  the  subscript  "c"  refers  to  the  center  of  the  earth,  u  is  the  angular  velocity  of  the  earth's 
rotation,  and  p*  is  a  vector  from  the  center  of  the  earth  to  the  radar  site.  We  have 


r  -4 

—  ~  10  ; 
c 

a 

—  -  1.6X10 


x  2  v  dv 

Af  2  dt 
c 


"fi  ~  10 

c 


-6 


-11 


SL.2  -  7  x  10 

c 


-11 


1  ~  to"4 


1  ~  10-10 


1 0  sec 


while  if  we  just  include  the  earth,  moon,  and  sun  in  the  gravitational  potential, 
,2  d|  /  iVAs  s  .  JVim  m  \  2 


Lf  dt 


(M  dr  M  dr  \ 

_ s  _ s  ,  _ m  _ m  \  T  < 

2  dt  2  dt  /  f 


(6-123) 


(6-124) 


where  rm  and  rg  are  the  distances  to  the  moon  and  sun,  respectively.  These  terms  do  not  ex¬ 
ceed  10"9  seconds.  Thus  we  can  use  Eq.  (6-121)  for  this  experiment. 
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VII.  THE  DOPPLER  EFFECT 


We  will  begin  by  outlining  the  theory  of  geometrical  optics  in  the  Special  Theory  of  Relativ¬ 
ity.58  We  look  for  solutions  of  Maxwell's  equations  where  the  quantities  describing  the  field,  say 
the  potentials,  are  of  the  form 

A  (x)  =  Re  [a  (x)  e^X^]  (7-1) 

where  p(x)  is  called  the  eikonal.  This  is  the  generalization  of  a  plane  wave  for  which 


A.(x)  =  Re  {[a.  exp[i(k.  x^  +  a)]}  =  Re  {fa.  exp[i(k  •  "r  -  wt  +  o?)]}  .  (7-2) 

We  want  solutions  of  the  form  of  Eq.(7-1)  which  are  approximately  plane  waves,  i.e.,  the 
amplitude  and  direction  of  the  wave  remain  practically  constant  over  distances  of  the  order  of  a 
wavelength.  Since  ip(x)  changes  by  2?r  in  one  wavelength,  ip(x)  will  be  a  large  number  for  a  macro¬ 
scopic  distance  and  furthermore,  derivatives  of  ip  will  be  much  bigger  in  size  than  derivatives 
of  a.(x).  Substituting  Eq.(7-1)  into  Maxwell's  equations,  and  making  use  of  the  assumed  large 
value  of  ip(x),  we  obtain  the  eikonal  equation: 


ij  dp  _  0 

0xl  axj 


(7-3) 


which  we  rewrite  in  terms  of  the  wave-vector: 

k  def_3£  kl  =  (k,co)  ;  7,..kV=0  .  (7-4) 

1  3xl  13 

A  surface  of  constant  p  is  called  a  wave-surface  and  the  wave-vector  is  normal  to  the  sur¬ 
face,  with  its  direction  defining  that  of  the  ray  (a  curve  whose  tangent  at  each  point  coincides  with 

the  direction  of  propagation  of  the  wave).  Equation  (7-3)  looks  like  the  Hamilton- Jacobi  equation 
59 

for  a  free  particle 

ij  _3S  +  m2  =  0  (7-5) 

ax1  ax3 


for  which  the  momentum  and  energy  of  the  particle  are  given  by 


pl  =  (?,  E)  ;  =  -m2  (7-6) 

so  that  a  ray  looks  like  the  trajectory  of  a  massless  particle.  The  analogy  can  be  made  even 
stronger  if  we  observe  that  Hamilton's  equations 

(7-7) 


(7-8) 

so  that  in  a  vacuum  where  |k|  =  w  the  speed  of  light  along  the  ray  is  unity  and  the  rays  are 
straight  lines. 


dP 

dt 


^rr  =  ~VE 


V  -  —  =  ^ 
dt  DrT 

ap 


have  as  analogues 


dk 

dt 


-rr  =  -Vco  ;  -rr  = 


dr 

dt 


dOJ 

dk 
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Now,  by  its  definition,  the  ray  is  perpendicular  to  the  wave  surface,  i.e.,  if  dxl  represent 
an  infinitesimal  vector  connecting  two  points  on  the  wave  surface  then 

0  =  dip  =  -^-r  dx*  .  (7-9) 

ax1 

But  the  ray  also  has  the  peculiar  property  of  lying  in  the  wave  surface,  i.e.,  if  we  consider  an 
infinitesimal  vector  along  the  ray  (ck1),  then  it  connects  points  on  the  wave  surface.  To  prove 
this,  note  that 

dip  =  — ek*  =  er — — — -r  =  0  (7-10) 

3xl  3x^ 

Let  us  now  look  at  an  emitter  and  receiver  of  light  and  consider  two  rays  joining  them/  ^  In 
Fig.  3  we  have  drawn  a  space-time  picture  of  the  process.  Note  that  ip  has  the  same  value  at 
each  point  along  the  ray  since  the  ray  lies  in  the  wave  surface.  We  now  define  the  "unit"  tangent 
vectors  to  the  world  lines  of  the  emitter  and  receiver 


dx, 


(1) 


dx 


(1)  ds(i)  '  U)  ds 


ds(D  =~’Jijdx(1i)dx(1i)  : 

u(ll)u(Vij  =  -1  =  u<2)u(Vij 


(2) 

(2) 

\\  =  -n;,dx/7.dx,J, 


ij  (2)  (2) 


and  note  that  we  can  then  compute  A  ip  in  two  ways: 


(It)  uAvAS/.v  =  A  Ip  -  (-^r)  uLAS/?v 

W'd)  (1)  (1)  '3x '  (2)  (2)  (2) 


(7-11) 

(7-12) 


(7-13) 


or 


k(ll)U(l)iAs(l)  "  k(2)u(2)iAs(2) 


(7-14) 


or 


As 

As 


2 


^Vtz) 


(7-15) 


1  3-63- 5<20] 


Fig.  3.  Space-time  diagram  for  Doppler  effect 
and  red  shift. 
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But  As  is,  by  definition,  just  the  proper-time  interval,  i.e.,  the  time  interval  measured  by  a 
clock  carried  by  the  observer,  and  Aip  is  just  2?r  times  the  number  of  waves  each  observer  de¬ 
tects  in  this  time  interval  and  is  the  same  for  both.  Thus  the  ratio  ot  the  frequencies  as  meas¬ 
ured  by  the  receiver  and  the  transmitter  is 


f2  7?/As2  ffl 
'•  a‘2 


<klui)(2) 

^Vd) 


(7-16) 


This  can  be  rewritten  in  a  more  familiar  form  if  we  note  that  the  ratio  in  Eq.  (7-16)  is 
written  in  terms  of  Lorentz-invariant  quantities  and  so  can  be  evaluated  in  any  coordinate  system. 
Choose  a  system  in  which  the  transmitter  is  at  rest.  Then 


dx(°L)  "  ° 


dx, 


(1) 


ds 


=  1 


dt 


(1) 


ds 


(1) 


kiU(l)  =  k4 


,  i  ,  dx(  2)  dt 

kiu<2)  =  kt  ir 


(k 


(2) 


v(2)  +  k4) 


1  —  V, 


2  ,l/2 

(2)1 


k.k1  =  0  =  | k|2  —  kf  . 


k4  =^44k  =  —to  =  —  |k| 


h  _  1  ~n  •  v(2) 

ft  [1  ~v(  2)1 1/2 


(7-17) 


The  above  procedure  may  seem  like  a  rather  long-winded  way  to  obtain  Eq.  (7-17)  but  it  has 
the  advantage  of  being  usable  almost  word  for  word  in  GRT.58  The  GRT  generalization  of 
Maxwell's  equations  in  free  space  also  admits  an  eikonal-type  solution,  only  now 


ij  3 »  3if<  _  Q 

ax1  ax3 


(7-18) 


which  we  can  again  rewrite  as 

k  dlf  .§£  k  =  o  .  (7-19) 

1  ax1  1  3 

Proceeding  exactly  as  before,  but  with  g..  replacing  77..,  we  again  end  up  with  the  ratio  of 
proper  time  frequencies 


^2  =  (gj.jkluJ)(2)  (7-2 

fl  (gijklu3)(1) 

Before  we  can  use  this  formula,  however,  we  need  the  relationship  between  and  ^^)* 

special  theory,  these  quantities  are  equal  and  the  rays  are  null  geodesics,  i.e.,  straight  lines 
whose  tangents  have  zero  length.  This  suggests  the  GRT  result  If  we  define  a  curve  by 
means  of 


dx^ 

d\ 


=  g 


ij  dip 

dxi 


(7-21) 


then  a  straightforward  computation  shows  that  because  of  Eq.  (  (-19) 
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(7-22) 


d^x*  r  i  dx^  dx^  _  n 

dxz  jk  dX  dX 

dx^  dx^  _  n 
^ij  dX  dX 


so  that  if  we  know  the  equation  of  the  null  geodesic  connecting  the  transmitter  to  the  receiver 
and  the  special  parameter  along  it,  we  can  relate  and  k^  by  means  of  Eq.  (7-21): 


£  ■  fa  £  4«>/fa  #  4. 


(7-23) 


Fortunately  for  us,  the  special  parameter  is  easy  to  find  when  the  metric  is  "static,”  i.e.,  there 
exists  a  coordinate  system  in  which  the  g. .  do  not  depend  on  time  and  in  which  g^  =  g^  =  g^  =  0 
(Ref.  62).  This  can  be  seen  from  the  variational  principle  in  Eq.  (6-10)  for  the  static  case 

J 


a  f  f  dx^  dx^  _dt_  dt_ 

(  J  l  do-  da  ^44  dcr  da 


da  =  0 

The  Euler- Lagrange  equation  for  6t  immediately  leads  to 


d  /  dt  ,  _  . 
da  ^44  da 


or 


da 


<-g44> 


dt 


(7-24) 


(7-25) 


since  changing  a  by  a  constant  factor  changes  nothing. 
We  can  then  finally  write 


fn 


1 


1 

dx* 

dX(2) 

/ 

1 

dxl 

dx<V 

gij  “g44 

dt 

dS<2). 

(2)  / 

gij  -g44 

dt 

dsd)J 

(1) 


(7-26) 


This  formula  obviously  reduces  to  the  special  relativistic  one  when  g. .  =  77. ..  Suppose  now  that 
both  observers  are  at  rest.  Then 


dx(1D  =  dx<2i) =  dxa) =  0 


ds 


(2) 

f. 


(2) 


(2) 


(1) 

(2) 


_g44(l) 

(2) 


1/2 


dt 


f 


1  .  ,-l/2 

g44  —  g. .  g44' 


■(= 


-g44 

g44(l)\ 
g44(  2) ' 


(2)  / 


1  ,  >-l/2 

S44  _  g  ^  ^44 


g44 


(1) 


1/2 


(7-27) 


which  is  just  the  usual  red-shift  formula.  GRT  thus  gives  us  a  combined  " Red- Doppler"  fre¬ 
quency  shift. 

We  now  turn  to  the  implications  of  Eq.  (7-26)  for  an  interplanetary  radar  experiment.  We 
send  out  a  signal  at  frequency  f^,  it  is  received  and  reflected  by  the  planet  at  frequency  f^,  and 
returns  to  the  earth  with  frequency  f^.  We  are  interested  in 


f3-'i  =  ri(77  -  ‘) 


and  we  want  f^  —  f^  to  one  cps  if  f^ 
can  write 


1 0 

10  cps,  so  we  need  the  ratio  up  to  one  part  in  10 


(7-28) 

“10  \i  • 

We 
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i  dx(3 ) 


£3  £2 

f2  fl 


gij(3>  -g44(3>  q0)  ds(3) 
_ 1  i 

gij(2)  -g44(2)  q<2)  ds(2) 


1  dx<2) 


gij( 2 )  -g44(2)  k<2)  ds 


(2) 


_>_  k»  ±LU 

gU">  -844(1)  dsd) 


i  def  / dxl\ 

q  =  van  2- 

Consider  first  f2/f^: 


3 

photon 


i  def  /  dx1  \ 

c  =  IdT/i 


2 

photon 


(7-29) 


ds(l) 

dx  A,  x  * 

,  1  (2) 

g44(l)  dt 

gi.l(2)K(2)  dt 

dS(2) 

dxA  v 

lr  1  (!) 

g44(2)  dt  J 

gij(l)k(l)  dt  J 

def 

=  Mi  1 12 


For  the  first  factor, 

ds 

dt 


dx^  dJ}  1/2 

5ij  dt  dt 


ax  I  r  a 

—  dr|  =  is44-sa^  y‘ ] 


Gt  /?, 1/2 
v  w  ] 

=  [1  +  ip  —  ( 1  —  ip)  v2]1//2  =  [1  +  ip  —  v2]  *  ^2  +  0 (0v2) 


where 


O^v2)  =  0(^-v2)  =  0(v4)  ~  10‘16 

(1  +#t)  [1  +*4  -v,2]1^2 

[  1  1  =  2  1/2 
1  (1  +</,2)  [1 

1  +  Ip 2  -  (1  -  Ip^)  k2  •  v2 


1  +  ip£  1  —  (1  -  2j/»2)  k2  •  v2 


1  +  -  (1  -  ip^)  k4  •  vi  1  +  ipJ[  1  ~  (1  ~  2^)  kd  •  vi 


+  0(10"20) 


So 


[1  +  * 


2,1/2  4 

—  v .  1  1  —  ( 1  — 


1  V1 


*1  [1 


7T  2.1/2  , 

+  ib->  —  1  1 


a  -2 »2>  k2  •  ?2 

■  (1  -  2i/,1 )  kd  •  Vj 


(7-30) 


(7-31) 


(7-32) 


It  is  now  convenient  to  shift  to  the  coordinate  system  used  earlier  where  the  ray  travels  in  the 
x-y  plane  and  the  x-axis  is  directed  along  the  line  between  the  earth  and  the  planet.  Then  from 
Eqs.  (6-61)  to  (6-62): 
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1  +  ^  -  vl 


1  +i^2_v2 


1  +  ^1  “  V1 


1/2  1  -  ^1  + 


4M 


2M\  .  ...  , 

— )  v~>  +  y(t 0)  v0 

r1  /  2x  |  2  2y 


1  +  “  v2 


1/2  1  -  v,  - 


2x"  [if  V2x  +  f(t2>  V2y] 


1  —  v  — 
lx 


2M  X  */+  \ 

-  V,  +  yd. )  v. 

r1  lx  1  ly 


-20 

+  0(10  ^U) 


1  +  ^1  “  V1 

1  +  “  vl 


2i 


1/2 


1  -  v 


2x 


1  —  v 


lx 
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Finally,  if  we  perform  a  similar  analysis  on  f3/f2, 
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(7-34) 


where  e^  and  are  unit  vectors  along  the  outgoing  and  returning  rays,  respectively.  The 
square  root  can  be  ignored  since 
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where  the  first  term  is  just  the  one  predicted  by  Special  Relativity. 
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If  we  calculate  A  2  and  A^,  they  are  equal  and  in  fact 
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(7-35) 


A,  =  A,-  =  -  ~  (6t)  .  +  0(— 6tv) 

12  23  dt  out  r 


(7-36) 


where  from  Eqs.  (6-105)  to  (6-109) 
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[Eq.  (6-105)] 
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which  near  superior  conjunction  becomes 
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If  we  treat  the  earth  and  the  planets  as  coplanar, 
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[Eq.  (6-108)] 


[Eq.  (6-109)1 


The  verification  of  Eq.(7-36)  is  straightforward  but  rather  tedious  so  we  will  just  outline  the 
steps  based  on  Eqs.  (6-60)  to  (6-62)  and  the  definition  of  in  Eq.  (7-32).  First  we  note  that 
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We  now  need 
that 


y( 1 1 )  and  yft^).  Since  the  y-values  are  the  same  at  r^  and  r^,  Eq.  (6-60)  tells  us 
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which  we  approximate  to 
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and  then  use  Eq.  (6-62)  to  obtain 
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(7-39) 


(7-40) 


Finally  we  must  express  v  as  a  dot  product,  just  as  we  did  with  v  .  From  Fig.  4,  we  see  that 

y  ^ 

the  x-axis  is  in  the  direction  of  e .  0  while  the  z-axis  is  in  the  direction  of  r  X  rp.  Since  the 

12  p  E 

x— ,  y-,  and  z-axes  form  a  right-handed  coordinate  system. 
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Fig.  4.  Geometry  of  special  coordinate  system  used 
in  discussion  of  Doppler  shift. 
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Vx  V  *  e12 


(7-42) 


We  can  now  eliminate  all  reference  to  our  special  choice  of  axes  and  verify  Eq.  (7-36). 
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Vm.  SUMMARY 


We  have  tried,  in  this  report,  to  discuss  all  those  aspects  of  the  General  Relativity  Theory 
which  can  be  tested  at  present  or  in  the  near  future  by  performing  interplanetary  radar  experi¬ 
ments  and  then  combining  the  radar  data  with  data  deduced  from  optical  observations  to  produce 
ephemerides  for  the  planets.  The  anomalous  time  delay  and  the  corrections  to  the  Newtonian 
equations  of  motion  will  certainly  be  easy  to  detect,  but  the  corrections  to  the  Doppler  shift  will 
not  be  so  easily  verifiable.^ 
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